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Direct  nuclear  demagnetization  was  used  to  coo!  high-density  solid  ''He  into  the  magnet- 
ically ordered  phases.  The  magnetic  state  of  the  ^He  during  precooling  and  demagnetization 
was  indicated  by  pressure  change  detected  by  a  capacitive  strain  gauge.  The  entropy  of 
the  samples  was  reduced  by  precooling  to  temperatures  between  0.55  and  3.99  mK  in  a 
field  of  2.5  T.  For  a  bcc  sample  with  a  molar  volume  of  20.10  cm^/mole  and  remaining 
entropy  of  approximately  0.01i21n2  two  magnetic  transitions  were  observed.  A  similar  ex- 
periment on  two  hep  samples  with  molar  volumes  of  19.65  cm^/mole  and  19.65  cm^/mole 
was  carried  out.  Pressure  data  at  various  entropies  indicated  different  behaviors  near  zero 
field:  a  decrease  in  pressure  occurred  for  low-entropy  samples,  with  an  increase  in  pressure 
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for  high-entropy  samples.  Results  from  the  high-temperature  expansion  of  the  multiple- 
exchange  and  Heisenberg  models  can  explain  only  the  pressure  data  with  remaining  en- 
tropies 0.72i21n2  or  higher.  We  conclude  that  the  pressure  behavior  observed  in  samples 
with  entropies  of  0.46,  0.23  and  0.027?ln2  are  the  results  of  magnetic  ordering  in  hep  solid 
3He. 
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CHAPTER  1 
INTRODUCTION 

Solid  ^He  becomes  magnetically  ordered  at  temperatures  <  1  mK  via  exchange  inter- 
actions made  possible  because  of  its  large  zero  point  motion.  Despite  its  simple  structure, 
■^He  has  rich  magnetic  properties.  The  body-centered-cubic  (bcc)  solid  in  a  magnetic  field 
has  three  magnetic  phases:  the  paramagnetic  phase  (PP),  the  low  field  phase  (LFP),  and 
the  high  field  phase  (HFP). 

While  most  of  the  magnetic  properties  of  the  bcc  solid  are  well  characterized,  there  exists 
little  information  about  the  properties  of  the  hexagonal-close-packed  (hep)  solid  which  forms 
at  high  pressure.  Theoretical  calculations  [1]  and  experimental  works  [2,  3,  4]  indicate  that 
the  solid  would  order  ferromagnetically  near  20  //K.  Such  a  low  temperature  makes  the 
study  of  the  magnetic  properties  in  the  hep  solid  difficult.  Despite  many  attempts  [2,  3,  4], 
no  magnetic  transition  in  the  hep  solid  was  observed  because  the  ordering  temperature 
could  not  be  reached  using  the  indirect  cooling  method. 

To  achieve  lower  temperatures,  ^He  samples  may  be  directly  demagnetized,  an  idea  that 
has  been  around  for  many  years  [5].  This  method  was  recently  employed  at  the  University 
of  Tokyo  [6]  on  both  bee  and  hep  solids.  Their  results  showed  that  the  bcc  solid  reached 
both  HFP  and  LFP  implying,  that  the  sample  was  cooled  below  28  fiK  (extrapolated  from 
susceptibility  data  of  high  density  bee  solid  ^He  [7]).  For  their  two  hep  samples,  no  clear 
signature  of  the  ordering  was  observed  because  of  heating  of  the  cell  near  the  end  of  the 
demagnetization . 


Using  the  sample  demagnetization  technique,  we  carried  out  pressure  measurements  for 
one  bcc  and  two  hep  samples.  The  pressure  changes  of  the  samples  during  the  precooling 
and  demagnetization  indicated  the  magnetic  state  of  the  sample.  With  an  improved  design 
of  the  experimental  cell,  we  have  reduced  unnecessary  heating  near  the  end  of  the  sample 
demagnetization  and  thus  eliminated  the  problem  encountered  by  the  Tokyo  group.  Also, 
we  have  made  the  first  clear  observation  of  the  pressure  signature  of  magnetic  ordering  in 
the  hep  solid  ^He. 

This  thesis  is  arranged  into  five  chapters.  In  chapter  one,  a  review  is  given  of  the 
current  understanding  of  the  nuclear  magnetic  properties  of  solid  ^He.  There  will  be  a 
discussion  of  the  existing  theoretical  models  and  a  summary  of  previous  experiments  on 
solid  ^He.  Chapter  two  is  devoted  to  the  experimental  arrangements  used  in  this  work.  A 
detailed  description  of  the  experimental  apparatus  and  data-taking  process  will  be  provided. 
Chapter  three  is  about  the  experimental  procedures  such  as  calibration  of  the  strain  gauges 
and  temperature  regulation.  The  experimental  results  and  their  contribution  toward  our 
understanding  of  nuclear  magnetism  in  solid  "^He  are  presented  in  chapter  four  together  with 
the  simulation  data  of  sample  demagnetization  using  the  multiple  exchange  and  Heisenberg 
models.  Chapter  five  is  a  summary  of  this  experiment  with  a  proposal  for  future  works. 

1.1    Pressure  and  Magnetic  Phase  Diagrams  of  Solid  "^He 

Because  of  weak  nuclear  interactions,  magnetic  ordering  of  nuclear  spins  in  most  solids 
exists  only  in  the  nanokelvin  temperature  range  or  lower.  In  ^He,  however,  ordering  takes 
place  at  a  few  mK  for  some  density  and  magnetic  field.  Such  a  high  ordering  temperature 
is  brought  about  by  particle  exchange  made  possible  because  ^He  has  a  large  zero-point 
motion.  The  wave  function  of  one  atom  overlaps  with  those  of  its  neighbors  leading  to  the 


exchange  of  their  positions  in  the  lattice.  Based  on  nuclear  magnetic  resonance  (NMR) 
data  [8],  the  exchange  rate  in  ^He  is  estimated  to  be  as  high  as  20  MHz  at  melting  pressure. 

The  pressure  versus  temperature  (F  -  T)  phase  diagram  of  ^He  is  shown  in  Fig.  1.1. 
As  a  result  of  the  large  zero- point  motion,  ^He  remains  liquid  at  low  pressure  even  at  T  =  0 
and  does  not  solidify  until  its  pressure  reaches  almost  3  MPa.  Between  3  and  10  MPa, 
solid  ^He  has  a  body-centered-cubic  (bcc)  structure,  again  as  a  result  of  the  large  zero- 
point  motion.  At  melting  pressure,  the  solid  orders  into  an  antiferromagnetic  phase  with 
a  Neel  temperature  Tj^j  of  0.934  mK  [9]  at  zero  magnetic  field.  As  the  solid  is  compressed 
at  high  pressure,  it  is  more  difficult  for  the  atoms  to  exchange  their  places  because  of  the 
hard-core  potential.  The  difficulty  of  the  atoms  to  exchange  at  high  density  is  analogous  to 
the  difficulty  of  moving  around  in  a  crowded  bus;  the  denser  the  crowd,  the  harder  it  is  for 
two  passengers  to  exchange  their  positions  [10].  The  exchange  frequency  becomes  smaller 
at  high  pressure,^  causing  a  rapid  decrease  in  the  ordering  temperature,  which  drops  below 
40  /iK  near  10  MPa  at  zero  field.  There  is  a  simple  empirical  equation  relating  the  ordering 
temperature  to  the  molar  volume  of  a  bcc  sample  [11]:         .-  '  v 

s 

Tn  =  v'^-'^\  (1.1) 

where,  v  the  molar  volume  of  the  solid,  ranges  from  24.14  to  22.46  cm^/mole,  corresponding 
to  the  solid  pressures  from  3.52  to  5.32  MPa  respectively.  '  , 

Above  10  MPa,  the  crystal  structure  of  solid  ^He  is  hexagonal-close- packed  (hep).  Be- 
cause of  the  lattice  symmetry  which  favors  three-particle  exchange,  the  hep  solid  is  expected 
to  order  into  a  ferromagnetic  phase  [1].  Previous  experimental  results  also  indicate  a  strong 

'The  exchange  frequency  also  decreases  for  processes  which  involve  more  thcin  four  particles  due  to  a 
longer  path  length. 
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Figure  1.1.  P-T  phase  diagram  of  ^He  at  zero  magnetic  field  where  F,  AF  and  PP  indicate 
the  ferromagnetic,  anti-ferromagnetic,  and  paramagnetic  phases  of  the  solid  respectively. 
A  and  B  show  the  superfluid  phases  of  the  liquid. 
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ferromagnetic  tendency  [2,  3,  4].  An  hep  sample  near  the  bcc-hcp  boundary  line  is  expected 
to  order  at  about  20  /xK  [3].  Such  a  low  ordering  temperature  makes  the  study  of  the  mag- 
netic properties  in  the  hep  solid  difficult. 

Despite  many  attempts  to  reach  the  ordered  phase  either  by  the  indirect  [2,  3,  4]  or 
direct  [6]  cooling  methods,  no  clear  evidence  of  magnetic  ordering  in  the  hep  solid  has  been 
observed.  The  main  goal  of  our  work,  therefore,  is  to  cool  an  hep  sample  into  its  ordered 
state,  identify  the  ordering  signature,  and  characterize  the  ordered  state. 

Unlike  the  hep  solid,  the  magnetic  properties  of  the  bee  solid  are  well  characterized, 
particularly  at  melting  density  {vmelt  =  24.22  cm^/mole).  There  are  many  advantages  to 
studying  a  solid  sample  on  the  melting  curve.  With  the  presence  of  the  liquid,  the  solid 
density  is  more  uniform.  The  high  ordering  temperature  (about  1  mK)  can  be  easily  reached 
by  nuclear  demagnetization  refrigerators.  Excellent  thermal  contact  between  the  solid  and 
the  refrigerator  can  be  provided  by  the  liquid.  A  low  pressure,  about  3  MPa  compared  to 
10  MPa  of  a  hep  sample,  yields  higher  sensitivity  in  pressure  measurements  because  the 
diaphragm  of  the  pressure  gauge  can  be  made  thinner.^  Experiments  on  higher-density  bcc 
solids  [11,  12]  show  that  many  of  the  magnetic  properties  can  be  scaled  to  those  at  the 
melting  density  using  a  simple  relation  such  as  Eq.  1.1. 

Despite  its  simple  lattice  structure,  ^He  has  surprisingly  rich  magnetic  properties.  Fig- 
ure 1.2,  taken  from  Ref.  [13],  shows  the  magnetic-field-versus-temperature  {B  —  T)  phase 
diagram  of  solid  ^He  at  its  melting  density.  There  exist  three  magnetic  phases:  the  para- 
magnetic phase  (PP),  the  low-field  phase  (LFP),  and  the  high-field  phase  (HFP).  In  the 
LFP,  the  nuclear  spins  have  an  up- two-down-two  (U2D2)  structure  [14].  It  is  made  up  of 

^The  pressure  change  in  the  sample  is  detected  by  the  deflection  of  a  metal  diaphragm.  A  thin  diaphragm 
deflects  more  than  a  thick  one  at  a  given  pressure.  See  section  2.1  for  discussion  on  the  construction  of  the 
pressure  gauge. 
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Figure  1.2.  B-T  phase  diagram  of  bcc  solid  ^He  at  melting  pressure. 
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ferromagnetic  planes  normal  to  one  of  the  three  principal  cubic  axes,  with  spins  pointing  in 
one  direction  in  two  planes  and  pointing  in  the  opposite  direction  in  the  next  two  planes. 
This  periodic  alternation  of  spin  direction  continues  throughout  the  entire  magnetic  domain 
of  the  crystal.  When  the  applied  field  is  larger  than  the  critical  field,  Bd,  the  solid  enters 
the  high-field  phase.  At  melting  density,  Bd  is  about  0.45  Tesla  and  decreases  rapidly 
with  increasing  solid  density.  In  1991,  Fukuyama  and  co-workers  [12]  measured  Bd  for 
molar  volumes  ranging  from  24.06  to  22.45  cm^/mole  and  found  that  at  zero  temperature, 
5ci(0)  =  yi^.sio.s  ^jjg  HFP,  the  spins  are  believed  to  order  in  a  canted  antiferromag- 
netic  structure  which  is  composed  of  two  simple-cubic  sublattices  with  spin  directions  tilted 
away  from  the  direction  of  the  applied  field.  The  angle  <p  between  the  spins  and  the  field  is 
equal  in  both  sublattices.  It  decreases  with  increasing  field  and  eventually  becomes  zero  at 
the  HFP-PP  transition  at  the  second  critical  field  5c2(0)  at  zero  temperature.  At  melting 
density,  Bc2(0)  is  expected  to  be  about  22  Tesla  [15].  The  existence  of  two  magnetically 
ordered  phases  in  the  bcc  solid  is  believed  to  result  from  the  competition  between  different 
exchange  processes.  A  detailed  explanation  can  be  found  in  the  next  section. 

1.2    Review  of  Exchange  Theory  of  ^He  Magnetism       ^.  ;  '. 

Currently,  there  is  no  existing  theoretical  work,  including  the  well-known  multi-particle 
exchange  theory  [1,  16],  that  can  accurately  describe  all  the  magnetic  properties  of  solid 
^He.  In  any  case,  a  closer  look  at  this  theory  is  worthwhile,  since  it  is  the  only  model 
which  can  explain  most  of  the  complicated  characteristics  of  magnetism  in  solid  ^He.  A 
simple  model  for  describing  exchange  of  particles  is  the  Heisenberg  nearest-neighbor  model 
(HNN)  which  takes  into  account  only  the  pairwise  exchanges  of  an  atom  and  its  nearest 
neighbors.-^  It  can  be  easily  shown  that  a  two-particle  exchange  process  causes  ^He  to  order 
^In  the  hep  solid,  three-particle  exchange  processes  also  lead  to  the  HNN  model. 
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antiferromagnetically  at  low  temperature.  The  proof  and  a  discussion  of  the  HNN  model 
is  given  below. 

1.2.1    Heisenberg  Nearest-Neighbor  Model 

^He  atoms  are  fermions  with  spin  1/2.  The  Pauli  principle  states  that  the  wave  function 
of  two  fermions  must  be  antisymmetric  when  the  particles  are  interchanged.  In  other 
words,  two  identical  fermions  cannot  occupy  the  same  quantum  state.  Mathematically,  this 
principle  can  be  expressed  as: 

r'  ^  '  .  V'i2  =  -V'2i '  ;, ;  .  •  •':  (1.2) 

where,  is  the  total  wave  function  of  the  two  fermions.  To  check  the  validity  of  this 
equation,  one  can  make  the  assumption  that  the  two  fermions  can  occupy  the  same  quantum 
state.  Equation  1.2  becomes  V'12  =  -V'12  =  0,  i.e.  the  amplitude  of  the  wave  function  is 
zero. 

The  wave  function,  V'12,  can  be  broken  down  into  two  parts: 

V'12  =  <?i»12Xl2,  (1.3) 

where  <j>i2  is  the  spatial  wave  function  and  X12  the  spin  state  of  the  two  fermions.  For  the 
ground  state,  (j)i2  must  be  symmetric;  that  is  =  -\-(t)2i-  For  two  spin-1/2  particles,  the 
total  spin  function  Xi2  can  be  either:  |  tit2>,  ^(|  tii2>  +|  iit2>),  ^(|  til2>  -Uit2>), 
or  I  4-1^,2  >.  Since  <^i2  is  symmetric  and  the  total  wave  function  ^12  =  <t>i2Xn  must  be  anti- 
symmetric, the  Pauli  principle  requires  that  X12  be  antisymmetric.  Out  of  the  four  choices 
only  -^{\  ti4'2>  -|  4-1 12>)  has  the  correct  symmetry.  This  spin  function  yields  a  total 
spin  of  zero  which  implies  antiferromagnetism.  It  can  be  shown  that  for  an  even  number 
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of  particle  exchange,  such  as  2,  4,  the  resulting  ordered  state  is  antiferromagnetic.  For 
odd  number  of  particle  exchange,  the  process  favors  a  ferromagnetic  state  (see  articles  by 
Roger  et  al.  [1]  and  Thouless  [16]  for  more  detail).  It  is  believed  that  the  competition  be- 
tween these  antiferromagnetic  (even)  and  ferromagnetic  (odd)  processes  brings  about  the 
existence  of  the  LFP. 

The  Hamiltonian  of  solid  ^He  is  straightforward, 

onon  "I"  H^xchange  "I"  H Zeeman  ipole- 

(1.4) 

The  energy  scale  for  each  term  is:  Hphonon  ~  10  K,  H Exchange  ~  1  niK  (at  melting  pressure), 
H Zeeman  ~  0.7  mK  (at  1  Tesla),  and  Hoipoie  ~  100  nK.  Because  of  the  large  separation 
between  these  values,  only  the  exchange  and  Zeeman  terms  dominate  the  thermodynamic 
properties  of  the  solid  for  the  temperature  range  between  10  fiK  to  100  mK.  Equation  1.4 
can  be  reduced  to 

H  —  H Exchange  "I"  Hz  eeman  ^  (1"^) 

and  the  partition  function  can  be  written  as 

Z  =  trace  {e~^^}  =  trace  {e-^('^£;^<^'"'"9=+^z«'"»'»)},  (1.6) 

where  /3  =  l/ksT  with  the  Boltzmann  constant.  For  the  exchange  of  two  nearest 
neighbors,  the  Hamiltonian  has  the  form  , 


^Exchange  —     "^Jnn  ^  .,  Sj  '  Sj,  ,  •   -■  (1*7) 

•<i  ' 


which  commutes  with  the  Zeeman  Hamiltonian, 


Hzeeman  =  -lh/2Y^Si-Bo,  (1.8) 

i  ' 

where  Si  is  the  spin  operator  in  units  of  k  at  the  ith  lattice  site,  7  is  the  gyromagnetic  ratio 
with  a  value  of  7/27r  =  32.43  MHz/Tesla,  and  Bq  is  the  applied  field.  Therefore  Eq.  1.6 
becomes 

Z  -  trace{e"^(^^"'"'"9'^e~^(^2«""'"^}.  (1.9) 

To  study  the  effect  of  the  exchange  interaction  in  ^He,  one  need  not  cool  the  sample 
into  its  ordered  state.  Even  at  high  temperature  (T  |«/nn|/^B))  the  impact  of  the 
interaction  on  the  solid  is  large  and  can  be  measured  through  different  thermodynamics 
quantities  such  as  pressure,  heat  capacity,  and  susceptibility.  Mathematically,  the  effect  at 
high  temperature  can  be  shown  by  first  writing  down  the  relation  between  the  free  energy 
and  the  partition  function, 

F=--^lnZ.  (1.10) 

Since  the  temperature  is  assumed  to  be  high,  the  exponential  terms  in  Z  can  be  expanded 
in  powers  of  (3.  Equation  1.10  then  becomes 


ln2  +  e2^-e3^  +  ---+^(7^5o/3)2(l  +  fcB^^/3+^  +  ...)  + 


(1.11) 

where  is  Avogadro's  number  and  62, 63  and  0^  etc.  are  related  to  the  exchange  energy. 
For  the  Heisenberg  nearest- neighbor  model,  62  =  12J^„,e3  =  12J^„,  and  0^  =  xJnnf^kB, 
where  x  =  number  of  nearest  neighbors,  which  is  8  for  the  bcc  and  12  for  the  hep  solid.  To 
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obtain  the  value  of  pressure  of  the  solid,  one  takes  the  negative  of  the  derivative  of  F  with 
respect  to  the  volume,  giving 


The  heat  capacity  can  be  calculated  by  taking  the  second  derivative  of  the  free  energy  with 
respect  to  the  temperature.  In  zero  field  and  at  high  temperature,  one  finds 


a  =  -Tg  =  ^(e,/?^  -  63/?^  +  •••).  (1.13) 


The  susceptibility,  given  by  the  second  derivative  of  the  free  energy  with  respect  to  the 
magnetic  field,  is 

where  B  =  6lj-a2/8kg,  and  C  =  fioNA  {■yh/2)^  /ksv  with  /io  the  permeability  of  free  space 
and  V  the  molar  volume.  Using  the  HNN  model  gives  =  AJnn/kB  and  B  =  4(J„„/A;b)'^ 
for  a  bcc  solid. 

1.2.2    Multiple-Particle  Exchange  Model 

The  simple  pairwise  exchange  interaction  of  the  HNN  model  alone  cannot  explain  many 
of  the  experimental  results.  For  example,  the  model  predicts  the  ordered  state  of  a  bcc 
solid  to  be  a  normal  antiferromagnet  in  which  the  atoms  at  the  vertices  have  their  spins  up 
while  the  center  atoms  have  their  spins  down.  However,  NMR  [14]  and  neutron  scattering 
results  [17]  indicate  that  the  ground  state  actually  has  an  up-two-down-two  configuration. 
Another  example  of  the  insufficiency  of  the  model  is  that  the  ordering  transition  is  predicted 
to  be  second  order  and  to  occur  at  2  mK.  However,  in  1974  Halperin  et  al.  [18]  found  the 
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transition  to  be  of  first  order,  occurring  at  about  1  mK.  In  1975,  McMahan  and  Wilkins  [19] 
pointed  out  the  possibility  that  two-particle  exchange  might  occur  less  frequently  than  other 
exchange  processes,  their  reasoning  being  that  the  steric  effect  for  two  particle  exchange 
is  larger  than  other  exchanges.  Using  the  crowded  bus  analogy  mentioned  earlier,  one  can 
see  that  it  requires  more  room  for  two  passengers  on  the  bus  to  exchange  their  places  while 
less  room  is  needed  if  three  or  four  persons  agree  to  switch  their  places  in  a  cyclic  fashion. 
The  general  form  of  the  multi-particle-exchange  Hamiltonian  is  [1,  16] 


where  n  is  the  number  of  fermions  involved  in  an  exchange  process  and  a  represents  a 
specific  process  in  the  n-particle  exchange.  For  example,  n  equals  two  for  a  nearest-neighbor 
exchange  process.  The  two- particle  process  also  includes  the  second,  third-nearest-neighbor 
exchanges  which  have  different  exchange  frequencies;  a  is  therefore  needed  to  distinguish 
each  of  these  cases.  (-1)^  is  1(-1)  for  permutation  involving  an  even  (odd)  number  of 
particles.  F„  is  the  permutation  operator  of  n  particles  and  Jna  is  the  exchange  frequency 
for  each  permutation. 

For  an  exchange  process  involving  more  than  two  particles  (n  >  2),  both  directions 
(clock-wise  and  counter-clock-wise)  of  the  process  must  be  taken  into  account.  The  ex- 
change Hamiltonian  becomes 


where  V2  =  P2  and  —  (P„  -|-  P'^)  for  n  >  2.  The  different  directions  are  represented 
by  F„  and  P~*.  The  permutation  operator  F„  can  be  expressed  in  terms  of  spin  operators. 


(1.15) 


(1.16) 
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Godfrin  and  OsherofF  [20]  give  these  expressions  for  exchange  interactions  of  two,  three, 
four,  five  and  six  particles.  For  two-particle  exchange,  V2  has  the  form 

^2  =      =  ^(1  +     ■      =  \  +  25.  •  S,-,  (1.17) 

where  the  Pauli  spin  matrices  (T,-  are  replaced  by  the  spin- 1/2  operators  Si-  For  the  three- 
spin  exchange, 

V':,  =  Piik  +  P^l^\  +  2{Si-S^  +  SrSk  +  Sk-Si).        ■■  (1.18) 

Note  that  Pijk  is  just  another  way  to  write  P3.  For  an  exchange  process  involving  four 
particles, 

n  =  PnM  +  P-,h  =  ^  +  E      •  '^'^  +  "^G'Jkl,  (1.19) 
where  the  summation  runs  over  the  six  distinct  pairs  (//,  u)  of  the  four  spins        k,  I)  and 

Gijki  =  {Si  ■  Sj){Sk  •  Si)  +  {S,  ■  Sk)iSi  ■  Si)  -  {Si  ■  Sk){Sj  ■  S,).  (1.20) 

If  one  explicitly  writes  out  every  term  of  the  summation  then 

Pijki  +  P-,li  =  ^  +  S,-Sj  +  Sj  -Sk  +  Sk-Si  +  Si-  S,  +  5, .  Sk  +  Sj  ■  Si  +  (1.21) 
4[(5,  •  S,){Sk  ■  Si)  +  {S,  .  Sk){Si  ■  5.)  -  (5.  •  5,) (5,  •  S,)]. 

For  five-spin  exchange, 

^5  =  -Pijfc/m  + -f^jfc^/^  =  -  +  -        5^  •  5^  +  2  Gf,^r,i-  (1.22) 

Ai<i/  /i<i/<r)<e 
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Finally,  for  the  exchange  of  six  spins 


'P'e  -  Pijklmn  +  Pijllmn  -  '^Yl        '  ^"'^  G f^^^ric  +  "^Sijklmn,  (1-23) 


where 


S^,k^mn  =  [{Si  ■  Sj){Sk  •  5,)(5„  •  Sn)  +  {Sj  ■  Sk){S,  ■  Sm)iSn  '  5,)]  (1.24) 
-[(5.-  •  Sj){Sk  •  Sm)iSl  ■  Sn)  +  {Sj  ■  Sk){Sl  ■  Sn){Sm  '  5.)  +        •  S,){Sm  "  5.)(5„  •  5,) 
+  {Sl  ■  S^){Sn  ■  S,){S,  ■  Sk)  +  {Sm  ■  Sn){Si  ■  Sk){Sj  ■  S,)  +  (5„  •  S,){S,  ■  S,){Sk  ■  Sm)] 
+  [{S,  ■  S,){S,  ■  Sn)iSk  ■  Sm)  +  (5,  •  Sm){Sk  '  Si){Sl  ■  5„)  +  (Sk  ■  Sn){Sl  ■  S,){Sm  "  5.)] 

+[{Si  ■  Sj)iSk  ■  Sn){S, .  Sm)  +  {Sj  ■  Sk){Si  ■  Si){Sm  ■  5„)  +  {Sk  ■  5,)(5„  •  5,)(5„  •  5.)] 
-[{Si-Si){Sj-S,n){Sk-Sn)]. 


Substituting  all  of  the  listed  cycles  back  into  Eq.  1.16,  one  gets 


H exchange     —     ^  ^ 


-T,J^a(^\  +  2Si-S,^  (1.25) 


+     E  J3Jl  +  2{Si-Sj  +  Sj-Sk  +  Sk-Si)) 

i<3<k  ^ 

-        E     '^4a  ^  Sfj,  •  Si,  +  AGijkl 

i<j<k<l  \         H<i^  y 

+     E   J^A\  +  \Y.s,-s.  +  2 

i<j<k<l<m  \  M<f  n<i/<ri<e  J 

~  E  "^6"  iTfi'^lE'^'"''^''"'"       E      GfJ-'^Vc  +  '^Sijktmn 


,  16  4 


Another  way  to  express  the  above  Hamiltonian  is  to  break  exchange  cycles  into  two- 
spin  terms,  four-spin  terms...,  and  then  sum  each  type  of  the  two-spin  terms.  For  example, 
in  Fig.  1.3  (which  is  redrawn  from  Ref.  [21]  with  the  next-nearest-neighbor  process  J„„„ 
added),  the  three-particle  exchange,  Jf ,  can  be  broken  into  two  nearest-neighbor  terms,  and 
one  next-nearest-neighbor  term.  The  four- particle  cycle  Kp  can  be  divided  into  four  nearest- 
neighbor  terms  or  into  one  next-nearest-neighbor  term  and  one  third-nearest-neighbor  term. 
Note  that  Jt  and  Kp  are  the  notations  used  by  Roger  et  al.  [1]  for  particular  exchange  cycles 
in  and  J4a,  respectively.  See  Refs.  [20,  22]  for  the  contributions  of  other  exchange  cycles. 
For  some  calculation,  it  is  more  convenient  to  express  the  Hamiltonian  this  way  [20,  23]. 
Using  the  decoupling  scheme  with  the  constant  terms  removed  for  clarity,  the  Hamiltonian 
becomes 

(1)  ^    ^  (2)  ^    ^  W  ^  ^ 

Hexchange     —     — 2Ji  ^  5,  •  Sj  —  2J2  ^  -^i  '  5j  2Jk  ^  Si  •  Sj  —  ■  •  ■  (1.26) 

i<j  i<i  i<j 

+  0(8^) +0(8^), 

where  the  four-spin  term  0(8"^)  has  the  form  (5,  •  Sj){Sk  ■  Si)  and  0(5^)  =  (5,  •  Sj){Sk  ■ 
Si){Sm  •  Sn)-  The  summation  's  taken  over  the  k*^  nearest  neighbor  and  Jk  has  the 
form  J2na  ^naJna  with  c„a  a  Constant.  For  example,  it  can  be  shown  that  for  the  exchange 
cycles  given  in  Fig.  1.3 

Jl    =   Jnn-QJt  +  SKp  (1.27) 

J2     =     Jnnn  —  ^Jt  +  Kp 
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Figure  1.3.  Schematic  drawing  of  important  exchange  processes  in  bcc  solid  ^He. 


In  1983,  Roger,  Delrieu,  and  Hetherington  [1]  published  a  review  paper  in  which  they 
considered  three  and  four  particle  exchange  effects  in  the  magnetism  of  bcc  solid  "^He.  Their 
goal  was  to  find  a  Hamiltonian  which  gave  results  consistent  with  experimental  data  and 
at  the  same  time  contained  the  smallest  number  of  exchange  processes.  They  found  that 
three-  (Jt)  and  planar  four-particle  (Kp)  exchanges  would  give  reasonably  good  agreement 
with  existing  experimental  data.  The  exchange  Hamiltonian,  Eq.  1.25,  is  reduced  to  only 
two  terms 

Hexchange     =  Jt         +  2{S,  ■  Sj  +  Sj  ■  Sk  +  Sk  ■  Si))  (1.28) 

i<3<k  ^ 
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i<j<k<l         \  ^^<l^ 


As  in  the  HNN  model,  this  Hamiltonian  was  used  to  calculate  the  free  energy  (see 
Eq.  1.10).  Various  thermodynamic  properties  at  high  temperature  were  determined  from 
the  expansion  of  the  free  energy.  For  the  bcc  solid,  Roger  et  al.  [1]  found 

1  o 

0^=      {Kp-2Jt),  (1.29) 
Kb 

e2  =  (24)2(j2  _  Ik^j^  +  lli^2)  (1  30) 

The  expression  for  0^  clearly  shows  the  competition  of  three-  (ferromagnetic)  and  four- 
particle  (antiferromagnetic)  exchanges.  The  values  of  Jt  and  Kp  are  obtained  by  fitting 
to  the  experimental  results  of  the  high-temperature  coefficient  62  and  the  mean  spin-wave 
velocity  Vspin  which  are  found  to  be  7  (mK)^  [24]  and  8  cm/sec  [25]  respectively.  The 
expression  for  Ujp,„  used  in  the  fitting  is  [1] 

Vspin    =    2(^^j  (2Ji  +  J2-h4J3  +  3A'p)'/'(-J2)'/M  J2+4J3-^j  (1.31) 
=   4f(3/f,-4J0^/^(4J,-K,)i/^(8J,-^)i^  ^ 

where  a  is  the  lattice  constant  and  J,  are  given  in  Eq.  1.27  with  J„„  =  J„„„  =  0.  For 
bcc  solid  at  the  melting  density,  Jt  has  a  value  of  -0.13  mK  and  Kp  of  -0.385  mK  yielding 
=  -2.25  mK  and  Tjv  =  1.2  mK  (compared  to  the  experimental  value  of  0.934  mK 
measured  by  Ni  et  al.  [9]). 
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Using  the  two  exchange  model,  Roger  et  al.  were  able  to  predict  the  U2D2  structure 
in  the  ordered  low-field  phase,  a  canted  structure  for  the  high-field  phase,  and  a  first-order 
transition  line  between  the  paramagnetic  phase  and  the  high-field  phase  in  the  low-field 
region.  Although  this  is  a  much  better  model  than  the  HNN  model  it  also  has  weaknesses. 
For  example,  it  predicts  the  critical  field  between  the  LFP  and  the  HFP  to  be  1.6  Tesla  at 
zero  temperature.  In  reality  Bc(0)  was  measured  to  be  0.452  Tesla  by  Xia  et  al.  [13].  To 
lower  the  critical  field,  Roger  et  al.  introduced  a  third  exchange  process  into  the  Hamilto- 
nian.  This  process  is  the  exchange  of  two  nearest  neighbors  J„„  used  in  the  HNN  model. 
Unfortunately,  this  process  creates  an  extra  ordered  phase  in  the  phase  diagram,  which  has 
not  been  observed  experimentally. 

In  1987,  Ceperley  and  Jacucci  [26]  calculated  the  exchange  frequencies  of  various  pro- 
cesses using  the  path-integral  Monte  Carlo  method.  They  found  that  even  six-particle 
exchange  processes  are  significant  in  bcc  "^He. 

In  the  same  year,  Roger  et  al.  [23]  published  the  results  of  their  high-temperature 
expansion  with  a  multiple-exchange  Hamiltonian.  They  included  in  the  calculation  two-, 
three-,  and  four-particle  exchange  processes  for  both  bcc  and  hep  solids  in  an  arbitrary 
field.  The  expansion  was  carried  to  fourth  order  in  /?  =  l/ArgT,  one  order  higher  than  the 
1983  calculation  [1].  Their  exchange  Hamiltonian  has  the  form 


where  Ji,  J2  are  the  combination  of  two-  and  three-particle  exchange  as  they  appear  in 
Eq.  1.27  and  K  is  the  exchange  frequency  involving  four  particles,  with  Yl^^^  calculated 


(1.32) 


over  the  nearest-neighbor  pairs  and 


over  the  next-nearest-neighbor  pairs.  For  a  bcc 
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solid,  the  expressions  for  the  interjictions  are: 


Ji    —    Jnn  —  6Jt 

J2     =     Jnnn  -  4  Jj         ,       ■  (1.33) 


K   =  Kp. 


Note  that  Ji  does  not  include  Kp;  therefore,  it  is  not  equal  to  J,-  which  is  defined  in  Eq.  1.27. 
The  exchange  processes  are  shown  in  Fig.  1.3. 

For  an  hep  solid,  Roger  et  al.  [23]  took  into  account  five  dominant  exchange  cycles: 
Jnni  J'nni  Jti  J'ti  ^-id  K'^^.  These  five  dominant  processes  together  with  three  less  important 
exchange  processes  (J*,  J^p,  and  K'p)  are  drawn  in  Fig.  1.4,  which  is  taken  from  Okamoto's 
Ph.D.  thesis  [27]  and  modified  slightly  for  clarity.  The  cross-hatched  atoms  are  shown 
smaller  in  size  to  give  a  three-dimensional  perspective.  As  in  the  bcc  solid,  the  same 
expression  for  the  exchange  Hamiltonian,  Eq.  1.32,  is  used  but  with 

Jl  Jnn  ~  Jt  ~  ^Jfi 

J2     =     J'nn-^J'u  '  (1-34) 

K  =  K^. 

Due  to  the  difference  in  the  crystal  structures,  J,  for  a  bcc  and  hep  solid  are  not  the  same. 
Note  also  that  J2  in  the  hep  case  is  for  the  nearest-neighbor  pairs  out  of  the  basal  plane 
instead  of  the  next-nearest-neighbor  pairs. 

Once  again,  the  free  energy  is  calculated  from  the  Hamiltonian  and  expanded  in  powers 
of  yS,  similar  to  Eq.  1.11.  Various  thermodynamic  quantities  can  then  be  obtained.  For 
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Figure  1.4.  Schematic  drawing  of  important  particles  exchange  processes  in  hep  solid  "^He 
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example,  by  taking  the  volume  derivative  of  the  free  energy,  up  to  the  first-order  term  of 
l/ksT  the  pressure  of  the  bcc  solid  is 


P^JBJ^P'^^^^]^...,  ,1.35) 


dv  2      ksdv  2 


and  that  for  the  hep  solid  is 


In  the  above  equations,  R  is  the  gas  constant,  p  the  polarization^  which  is  given  by  p  = 
tanh(^7Bo/2A;Br),  and  Oyj  the  Weiss  temperature 


4Ji    +3J2    +18Kp   bcc  solid 
kaOu,  =  {  (1.37) 

3Ji    -I-3J2    +QK'sg    hep  solid. 


Since  formulae  1.35  and  1.36  are  valid  for  high  magnetic  field  at  high  temperature,  unlike 
Eq.  1.12,  which  requires  hjEo/ksT  <C  1,  the  volume  derivatives  of  6^,  Kp,  and  K'^^  can  be 
calculated  from  the  high-T  pressure  data  of  a  ^He  sample  (see  Sections  4.1.1  and  4.2.1  for 
more  discussion).  Also,  in  Section  4.2.3  we  have  presented  detailed  description  of  how  this 
model  is  used  to  simulate  a  sample  demagnetization. 


^For  p  =  0  hcJf  of  the  toted  number  of  spins  cire  pointing  in  the  same  direction  with  the  applied  field  and 
the  other  half  are  pointing  in  the  opposite  direction.  For  p  =  1  all  of  the  spins  are  pointing  in  the  direction 
of  the  applied  field. 
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1.3    Previous  Experiments  on  Solid  ^He 

1.3.1    bcc  Solid 

Because  of  the  simple  nature  of  ^He,  it  was  thought  that  the  ordinary  HNN  model  would 
explain  all  the  magnetic  properties  of  solid  ^He.  Early  high-temperature  experimental  data 
seemed  to  support  this  model  quite  well.  However,  in  1971,  Kirk  and  Adams  [28]  made 
pressure  measurements  of  various  samples  at  different  fields  and  found  that  there  was  a  large 
discrepancy  between  the  theoretical  and  experimental  values  of  the  Weiss  temperature  0^- 
This  inconsistency  yielded  the  first  sign  of  problems  with  the  simple  HNN  model.  Below 
are  brief  descriptions  of  important  experiments  in  the  study  of  magnetism  in  solid  ^He. 

In  1974,  Halperin  et  al.  [18]  discovered  the  magnetic  ordering  transition  in  the  bcc  solid 
at  melting  density.  At  the  transition  temperature  of  about  1  mK,  they  found  an  abrupt 
decrease  in  entropy  indicating  that  the  transition  is  first  order.  In  the  recent  measurement 
by  Ni  et  al.  [29],  the  entropies  just  above  and  below  the  transition  temperature  are  found 
to  be  0.63i21n2  and  0.19i21n2,  respectively.  The  discovery  by  Halperin  et  al.  contradicts 
the  prediction  by  the  HNN  model  of  a  second-order  transition  at  2  mK.  The  observation 
sparked  much  interest  and  many  experiments  soon  followed.  Realizing  the  insufficiency 
of  the  simple  HNN  model,  theorists  have  sought  better  models  which  can  explain  the  rich 
magnetic  behavior  of  solid  ^He.  One  of  these  models  is  the  two-exchange  {Jt  and  Kp)  model 
mentioned  in  the  last  section. 

In  1977,  Kummer  et  al.  [30]  carried  out  an  experiment  which  was  similar  to  the  one 
performed  by  Halperin  et  al.  but  in  a  finite  field  up  to  1.2  T.  They  found  that  the  ordering 
transition  decreased  with  increasing  field  up  to  a  field  of  0.4  T.  With  higher  field  the 
behavior  was  reversed,  an  increase  in  field  raised  the  ordering  temperature,  causing  a  kink 
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in  the  phase  diagram.  They  correctly  identified  the  kink  as  an  indication  of  a  new  ordering 
phase  which  is  known  today  as  the  HFP  (see  Fig.  1.2). 

The  transitions  to  the  LFP  from  HFP  and  PP  are  accepted  to  be  first  order.  However, 
there  were  disputes  about  the  order  of  the  PP-HFP  transition.  An  NMR  experiment  in  1981 
by  Osheroff  [31]  in  fields  up  to  0.6  T  indicated  that  the  transition  was  first  order.  Similar 
conclusions  were  also  drawn  from  two  other  works  [32,  33].  In  other  experiments  [34,  35],  the 
results  indicated  second-order  behavior.  This  disagreement  was  resolved  in  1993  by  Xia  et 
al.  [13].  In  this  work,  they  found  that  the  transition  was  first-order  for  B  <  0.65  Tesla, 
indicated  by  the  discontinuity  of  the  entropy  across  the  transition.  In  higher  fields,  the 
entropy  changed  continuously  through  the  phase  transition,  so  here  the  transition  is  second 
order. 

1.3.2    hcD  Solid  

While  many  significant  works  have  been  done  on  the  bcc  solid  at  low  temperatures,  there 
has  not  been  much  work  on  the  hep  solid  because  the  magnetic  ordering  of  the  hep  solid 
occurs  at  very  low  temperatures.  An  hep  solid  near  the  bcc-hcp  boundary  is  expected  to 
order  at  about  20  /iK  or  less.  The  first  attempt  to  cool  an  hep  sample  into  its  ordered  state 
was  made  in  1985  by  Takano  et  al.  [2].  Using  PrNis  and  copper  nuclear  demagnetization 
stages  they  cooled  their  samples  to  43  //K,  but  no  ordering  signal  was  observed.  However, 
by  measuring  the  AC  susceptibility  of  samples  with  molar  volumes  ranging  from  18.1  to 
19.4  cm^/mole  and  fitting  the  data  to  a  Curie-Weiss  form,  x  =  C/{T  -  B^j),  they  found  a 
positive  Weiss  temperature  O^j  for  all  samples.  The  result  is  an  indication  of  a  ferromagnetic 
tendency  in  the  hep  solid.  The  fitting  form  is  similar  to  the  Eq.  1.14  but  with  the  third 
term  in  the  denominator  set  to  zero  in  the  high-temperature  approximation. 
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A  similar  experiment  was  carried  out  by  Yano  et  al.  [3]  in  1990.  In  this  work,  the  samples 
were  cooled  to  38  /iK  by  double  nuclear  demagnetization  (PrNis  and  Cu  in  series).  Again,  all 
Weiss  temperatures  were  found  to  have  a  positive  value,  but  there  was  no  transition  found 
in  any  of  the  samples.  However,  the  molar  volume  dependence  of  6^  was  experimentally 
determined  to  be  v^^,  which  was  in  good  agreement  with  high  temperature  NMR  data  [8]. 

In  1989,  Miura  et  al.  [4]  made  the  first  pressure  measurements  on  an  hep  sample  down 
to  0.93  mK  in  magnetic  fields  up  to  0.2  T.  Fitting  the  results  to  the  high  temperature 
expansion  (Eq.  1.12)  and  using  only  de2/dv  and  dOu,/dv  terms,  they  found  d6^/dv  to  have 
a  positive  value,  another  indication  of  a  ferromagnetic  tendency  in  the  hep  solid. 

Even  with  the  best  cryostat  available,  current  technology  does  not  allow  high-density 
solid  ■'He  to  be  indirectly  cooled  into  its  ordered  phase  by  a  nuclear  demagnetization  stage. 
In  addition  to  the  low  ordering  temperature,  a  long  thermal  time  constant,  which  is  pro- 
portional to  1/T  in  zero  field  [35],  makes  it  essentially  impossible  to  reach  the  ordering 
transition.  Fortunately,  direct  demagnetization  of  solid  ^He  nuclei  (an  idea  that  has  ex- 
isted for  many  years  [5])  can  bypass  these  problems.  Recently,  this  method  was  used  at  the 
University  of  Tokyo  by  Okamoto  et  al.  [6]  in  an  attempt  to  cool  ^He  solids,  both  bcc  and 
hep,  into  their  ordered  phases.  This  group  successfully  cooled  a  bcc  sample  with  a  volume 
of  19.99  cm^/mole  into  its  HFP  and  LFP  phases.  However,  for  their  two  hep  samples,  with 
t;  =  19.66  and  19.07  cm^/moie,  respectively,  they  found  no  clear  ordering  transition. 

One  of  the  main  problems  with  the  Tokyo  experiment  was  heating  of  the  samples  near 
the  end  of  the  demagnetization,  attributed  to  an  indium  0-ring  which  became  supercon- 
ducting near  zero  field.  To  eliminate  the  problem,  we  have  designed  a  better  experimental 
cell  in  which  there  is  no  superconducting  0-ring.  The  description  of  the  cell  is  presented 
in  Chapter  2  along  with  details  of  the  experimental  setup  and  techniques. 


CHAPTER  2 
EXPERIMENTAL  APPARATUS 

The  objective  of  this  experiment  was  to  cool  hep  ^He  to  its  magnetically  ordered  state 
and  to  determine  the  ordering  by  measuring  the  thermodynamic  pressure.  In  order  to 
cool  an  hep  sample  into  its  magnetically  ordered  state,  the  sample  itself  was  used  as  a 
second  demagnetization  stage.  Our  hep  data  indicate  that  a  sample  temperature  of  20  fiK 
or  less  can  be  achieved  using  this  direct  cooling  method.  The  sample  in  a  field  of  2.5  T, 
along  with  a  copper  demagnetization  stage  in  a  field  of  7  T,  was  preeooled  to  11  mK  by  a 
dilution  refrigerator.  Then,  the  copper  stage  was  demagnetized  to  lower  the  temperature. 
Upon  reaching  the  desired  temperature,  the  sample  was  demagnetized  to  zero  field.  During 
the  sample  precooling  process,  the  temperature  of  the  copper  stage  was  measured  by  a 
^He  melting  pressure  thermometer  while  the  conditions  of  the  dilution  refrigerator  were 
monitored  by  various  resistance  thermometers. 

The  pressure  change  in  the  solid  was  detected  by  measuring  the  change  in  the  capac- 
itance of  a  strain  gauge.  Since  the  pressure  change  was  small  (on  order  of  10  Pa)  and 
the  measurements  were  made  at  microkelvin  temperatures  with  a  changing  magnetic  field, 
there  were  many  restrictions  in  the  designs  of  the  pressure  gauge  and  the  cell.  Nevertheless, 
with  a  new  design  of  the  cell,  we  have  avoided  significant  heating  encountered  by  the  Tokyo 
group  [6]  near  the  end  of  the  sample  demagnetization  resulting  in  a  better  understanding 
of  the  pressure  data  and  a  clear  observation  of  the  magnetic  ordering  in  hep  solid. 
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In  this  chapter,  the  experimental  apparatus  used  in  this  work  is  covered  in  detail.  The 
first  part  describes  the  experimental  cell  which  was  used  to  hold  the  sample  and  to  detect 
the  pressure  change.  Two  different  experimental  cells  were  made  for  this  work,  the  second 
cell  being  an  improved  version  of  the  first.  The  main  goals  for  the  second  design  were  to 
minimize  the  magnetic-field-dependence  and  maximize  the  thermal  conductance  of  the  cell. 
The  first  section  deals  exclusively  with  the  second  cell  in  which  magnetic  ordering  of  the 
bcc  and  hep  solids  was  observed.  The  second  section  describes  the  electronic  set-up  of  the 
bridge  used  to  detect  capacitance  changes  in  the  strain  gauge  of  the  cell.  The  third  section 
details  the  construction  of  a  2.5  T  home-made  magnet  which  was  used  to  demagnetize  the 
^He  sample.  The  fourth  section  describes  the  cryostat  and  the  fifth  section  details  the 
thermometric  methods  employed. 

2.1    Experimental  Cell 

A  schematic  drawing  of  the  experimental  cell  is  shown  in  Fig.  2.1.  The  pressure  change 
in  the  sample  was  detected  by  a  capacitive  pressure  transducer  of  simple  design.  One  wall 
of  the  sample  cell  is  a  thin,  flexible  diaphragm  which  moves  up  and  down  in  response  to 
increases  and  decreases  in  the  pressure  of  the  sample.  A  coin-silver  disk  is  attached  to 
the  diaphragm  and  serves  as  a  movable  capacitor  plate  while  another  disk,  held  close  to 
the  movable  plate,  is  used  as  a  fixed  plate.  As  the  pressure  changes,  the  gap  between  the 
two  plates  varies,  leading  to  changes  in  capacitance.  The  magnitude  of  the  change  can  be 
detected  by  a  capacitance  bridge,  the  design  and  operation  of  which  is  described  in  the 
next  section. 

Taking  pressure  measurements  during  the  precooling  and  nuclear  demagnetization  of 
high  density  solid  ^He  requires  two  main  considerations  in  the  design  of  the  cell:  the 
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Figure  2.1.  Schematic  drawing  of  the  experimental  cell. 
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maximization  of  the  thermal  conductance  between  the  sample,  the  cell  and  the  copper 
nuclear  stage  and  the  reduction  of  eddy-current  heating  which  exists  during  the  sample 
demagnetization.  While  the  first  condition  calls  for  a  large  area  in  the  design  of  the  cell, 
the  second  requires  a  small  cross-sectional  area.  Both  of  the  constraints  are  met  in  the 
design  of  this  cell,  as  described  in  the  following  text. 

The  upper  part  of  the  cell,  made  of  coin  silver  (90%  Ag  and  10%  Cu),  is  a  Straty- 
Adams  strain  gauge  which  has  a  fixed-edge  diaphragm.^  The  diaphragm  of  the  gauge  has  a 
diameter  of  1.22  cm  (0.48")  and  thickness  1.6  mm  (0.063")  chosen  so  that  the  strain  gauge 
can  withstand  high  pressure  and  at  the  same  time  yield  high  sensitivity.  The  formula  used 
for  this  determination  can  be  found  in  a  review  article  by  Adams[36]  and  is  given  by 


where  t  is  the  thickness,  a  the  radius  of  the  diaphragm,  the  maximum  operating  pressure 
and  Sm  the  maximum  stress.  For  this  cell,  Sm  was  taken  to  be  150  MPa  (21000  psi)  and 
Pm  was  15  MPa  (2100  psi). 

Learning  from  the  first  gauge  that  the  mylar  spacer  and  stainless  steel  screws  had  caused 
a  large  heat  load  and  magnetic-field-dependence,  the  second  gauge  was  designed  without 
them.  The  screws  were  no  longer  needed  because  the  holder  for  the  fixed  plate  had  threads 
on  its  inner  wall,  similar  to  the  cap  of  a  jar,  while  the  other  holder  had  similar  threads  on 
its  outer  wall.  These  were  machined  so  that  when  they  were  tightly  screwed  together  there 
was  a  0.05  mm  (0.002")  gap  between  the  two  plates  which  eliminated  the  need  for  a  spacer. 
The  gauge  was  thermally  linked  to  the  rest  of  the  cell  by  two  annealed  silver  wires  1.3  mm 
'The  diaphragm  is  machined  into  the  sample  cell  so  that  its  edge  is  held  fixed  by  the  waD  of  the  cell. 
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(0.05")  in  diameter  which  were  welded  to  the  holder  of  the  movable  plate  at  one  end  and 
to  small  silver  blocks  at  the  other.  The  blocks,  not  shown  in  the  figure,  were  tightly  fixed 
to  the  base  of  the  cell  by  stainless  steel  screws.  To  reduce  eddy  current  heating  during  the 
sample  demagnetization,  the  holder  of  the  fixed  plate  and  the  two  plates  had  slots  cut  into 
them.  The  use  of  the  annealed  silver  wires  together  with  the  presence  of  slots  at  various 
positions  on  the  cell  enable  the  two  constraints  of  maximizing  thermal  conductance  and 
minimizing  eddy-current  heating  to  be  satisfied. 

Electrical  wires  were  soldered  to  the  plates  using  EutecRod  #157  soft  solder.^  To 
protect  the  solder  joints,  the  wires  were  placed  inside  short  pieces  of  teflon  tube  and  a 
small  amount  of  Stycast  1266  epoxy"^  was  applied  to  hold  the  tube  in  place.  The  wires  were 
single  filament  superconducting  wires,  the  same  type  used  for  winding  the  1  T  magnet.'*  The 
two  capacitor  plates  were  glued  to  their  holders  with  Stycast  1266  epoxy  using  thin  porous 
Kimwipe  paper  wetted  with  Stycast  to  provide  electrical  insolation.  After  the  Stycast  had 
dried,  the  plates  were  re-machined  together  with  their  holders  to  ensure  that  they  were 
parallel  to  one  another. 

The  lower  part  of  the  cell,  made  of  silver  of  99.999%  purity,  was  a  cylinder  with  diameter 
of  1.3  cm  (0.5")  and  length  of  7.4  cm  (2.9").  A  shallow  flat-bottomed  hole  was  bored  in 
the  top  end  with  seven  shallow  holes  added  at  the  bottom  of  the  counter-bore  to  provide 
a  stronger  bond  between  the  silver  powder  and  the  cylinder.  Two  slots  of  dimensions 
0.05  cm  x  4.8  cm  (0.02"  x  1.9")  were  cut  into  the  body  of  the  cylinder  and  at  right  angles 
to  each  other  in  order  to  reduce  eddy  current  heating.  The  bottom  end  of  the  cylinder  was 
welded  to  a  silver  plate  of  thickness  of  0.31  cm  (0.125")  with  11  clearance  holes  for  screws 

^Eutetic  Corporation,  4040  172nd  St.,  Flushing,  NY  11358,  phone:  (800)  323-4845. 
^Deanco,  7130  University  Blvd.,  Winter  Park,  FL  32792,  phone:  (407)  671-5100. 

^To  ensure  that  they  would  remciin  superconducting  in  high  fields  and  thus  generate  less  heat  during  the 
pressure  meeisurements. 
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which  were  used  to  clamp  the  cell  to  the  copper  nuclear  demagnetization  stage.  Slots  were 
cut  between  these  holes  to  maximize  the  contact  area  between  the  cell  and  copper  stage  by 
allowing  the  silver  plate  to  flex.  Without  the  slots,  the  rigid  base  of  the  cell  could  have  as 
few  as  three  contact  points  to  the  copper  stage. 

The  lower  part  of  the  cell  was  annealed  for  12  hours  at  750°C  (1380°F)  with  a  pressure 
of  10~^  Torr  of  O2  gas.  To  reduce  the  contamination  from  oil  mist  and  other  foreign 
materials  inside  the  oven,  greened  stainless  steel  sheets  were  used  to  shield  the  silver.  After 
annealing,  the  silver  was  left  inside  the  oven  to  cool  to  room  temperature  over  night.  A 
residual  resistance  ratio  (RRR)  of  about  400  was  achieved  using  this  procedure.  A  copper- 
nickel  (Cu-Ni)  fill  line  with  a  diameter  of  5  mm  (0.02")  was  soft  soldered  to  the  body  of 
the  cell  1.59  cm  (0.63")  below  the  sinter  area  using  EutecRod  #157  solder. 

Silver  sinter,^  of  700  A  particles  pre-heated  at  295°C,  was  then  packed  in  the  cell  to 
provide  good  thermal  contact  between  the  bulk  ^He  and  the  rest  of  the  cell.  During  the 
packing  process,  a  packing  pressure  of  6.9  MPa  (1000  psi)  was  applied  to  a  small  amount 
of  sinter  which  was  added  at  each  step  to  ensure  the  same  density  throughout  the  cell.  The 
final  packing  ratio  of  50%  yielded  a  surface  area  of  approximately  1  m'^,  for  the  0.5  g  of 
sinter  powder. 

The  heating  problem  encountered  by  Okamoto  et  al.  [6]  near  the  end  of  their  sample 
demagnetization  was  caused  by  flux  jumps  in  a  superconducting  indium  0-ring.  The  best 
way  to  solve  this  problem  is  to  completely  eliminate  the  indium  ring  from  the  structure  of 
the  cell.  This  was  achieved  in  the  design  of  our  cell  by  replax;ing  the  indium  with  an  epoxy 
seal. 

^ULVAC  Technologies  Inc.,  6  Riverside  Dr.,  Andover,  MA  01810,  phone:  (508)  68&-7550,  fax:  (508) 
689-6300. 
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The  strain  gauge  and  the  body  of  the  cell  were  held  together  and  vacuum  sealed  by 
Stycast  2850FT  epoxy  prepared  with  low  viscosity  catalyst  24  LV.  This  catalyst  was  chosen 
because  it  has  better  adhesion  and  higher  resistance  to  thermal  shock  than  catalyst  9  or 
11.  The  epoxy  was  left  to  cure  for  about  12  hours  at  room  temperature,  as  recommended 
by  the  manufacturer.  The  cell  was  leak-checked  at  liquid  Nitrogen  (LN2)  temperature  with 
a  pressure  up  to  15MPa  (2200  psi),  safely  above  the  pressure  14MPa  (2000  psi)  at  which 
the  low  density  hep  solid  forms.  -      *  - 

2.2    Pressure-Measurement  Method        »  f "    '      ,  » 

The  capjicitance  change  in  the  strain  gauge  was  detected  using  a  bridge  configuration  in 
which  the  unknown  capacitance  of  the  strain  gauge  was  compared  to  a  reference  capacitor. 
The  reference  capacitor,  made  from  beryllium-copper  (BeCu),  was  designed  similar  to  the 
experimental  strain  gauge.  The  plates  of  the  capacitor  were  glued  to  their  holders  using 
Stycast  epoxy,  which  were  separated  by  a  mylar  spacer  with  a  thickness  of  0.05  mm  (0.002") 
and  bolted  together  to  achieve  mechanical  stability. 

For  high-precision  measurements,  the  reference  capacitor  must  be  stable,  bCjC  <  10~^, 
a  condition  which  is  almost  impossible  to  meet  at  room  temperature.^  Fortunately,  stability 
can  be  achieved  by  placing  the  capacitor  on  the  mixing  chamber  of  the  dilution  refrigerator. 
In  the  millikelvin  regime  the  thermal  expansion  coefficient  of  BeCu  is  small  so  that  any 
change  in  the  temperature  of  the  reference  capacitor  does  not  create  a  detectable  change 
in  its  capacitance. 

Three  coaxial  lines  were  used  in  this  electronic  arrangement.  The  strain  gauge  and 
reference  capacitor  each  had  their  own  isolated  line,  and  the  third  line  was  common  to 

*  Because  the  large  thermal  expansion  of  materials  at  high  temperatures  means  that  the  temperature  of 
the  capacitor  would  have  to  be  kept  extremely  stable. 
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both.  Below  the  1  K  pot  the  coaxial  cables  were  constructed  by  inserting  niobium-titanium 
(NbTi)  superconducting  wires  inside  CuNi  capillaries.  Before  insertion,  the  wires  were 
coated  with  vacuum  grease  to  keep  them  rigid  inside  the  tubes.  The  cables  were  thermally 
anchored  to  various  parts  of  the  refrigerator,  such  as  the  1  K  pot  and  mixing  chamber.  To 
improve  the  heat  sinking  of  the  coaxes,  the  thermal  contact  area  of  the  tubes  was  increased 
by  wrapping  the  tubes  around  Cu  formers  and  soft  soldering  them  together.  At  the  top  of 
the  cryostat  the  lines  were  terminated  inside  a  vacuum-tight  brass  header. 

A  schematic  of  the  bridge  is  shown  in  Fig.  2.2.  The  Eaton  1011 A  ratio  standard'^ 
and  capacitors  were  isolated  from  ground  by  a  1:1  Eaton  ST-200A  isolation  transformer. 
Referring  to  the  figure,  Cp  denotes  the  unknown  capacitance  of  the  pressure  gauge  of  the 
experimental  cell,  and  Cfle/  the  value  of  the  reference  capacitor.  A  variable  resistor  was 
used  to  balance  the  resistive  part  of  the  signal.  The  resistor  had  a  maximum  value  of  100  Q. 
and  was  placed  in  an  aluminum  box  to  shield  rf  noise. 

A  personal  computer  (316S  CompuAdd*)  was  used  to  control  the  output  of  the  HP  3325A 
function  generator^  and  to  read  and  store  data.  A  low  noise  pre-amplifier  (Ithaco  1201^°) 
was  added  to  improve  the  signal-to-noise  ratio.  The  output  signal,  Vb,  was  amplified  by 
10"*  before  being  fed  to  the  PAR  124A  lock-in  amplifier. In  addition,  the  high  pass  and 
low  pass  rolloffs  of  the  pre-amplifier  were  set  at  100  Hz  and  1000  Hz,  respectively,  to  nar- 
row the  bandwidth.  The  pre-amplifier  was  connected  to  the  coaxial  header  on  top  of  the 
cryostat  via  a  short  coaxial  cable  and  enclosed  in  a  polystyrene  box  to  protect  it  from  room 

^Eaton  Corporation,  Electronic  Instrumentation  Division,  5340  Alia  Road,  Los  Angeles,  CA  90066, 
phone:  (213)  822-3061. 

^CompuAdd,  12303-G  Technology  Blvd,  P.O.  Box  200777,  Austin,  TX  78727-0777,  phone:  (800)  999- 
9901. 

^Hewlett-Packard,  P.  O.  Box  105005,  Atlanta,  GA  30348,  phone:  (800)  452-4844,  fax:  (404)  980-7256. 
'"Ithaco,  Inc.,  735  West  Clinton  Street,  P.  O.  Box  6437,  Ithaca,  NY  14851-6437,  phone:  (607)  272-7640, 
fax:  (607)  272-4550. 

"EG&G  Princeton  Applied  Research,  P.  O.  Box  2565,  Princeton,  NJ  08645-2565,  phone:  (512)  272-9621 
or  (512)  272-9649. 
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Figure  2.2.  Electronic  circuit  of  the  capacitance  bridge. 
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temperature  fluctuations.  The  in-phase  signal  of  the  lock-in  amplifier  was  recorded  by  a 
Linseis  chart  recorder^'^  and  a  Hewlett-Packard  HP  3478A  multimeter  digitized  the  output 
signal  for  storing  in  the  computer. 

Except  for  the  computer,  everything  shown  in  Fig.  2.2  was  located  inside  the  shielded 
room.  Communication  between  the  computer  and  the  shielded-room  electronics  was  made 
through  two  HP  37204  HP-IB  extenders  which  were  placed  at  the  two  ends  of  the  commu- 
nication line  and  linked  together  by  a  fiber-optic  cable  in  order  to  eliminate  rf  noise.  The 
equipment  and  the  computer  were  connected  to  the  extenders  with  GPIB  cables. 

For  an  ideal  circuit,  the  input  from  the  function  generator  and  the  output  voltage  of 

the  bridge  are  related  by:  ^ 

Vb  _  {R-l)Cp  +  RCR,j^  •  '  ^2.2) 

Vi  Cp  +  Cuef 

where  R  is  the  value  of  the  ratio  standard.  When  the  bridge  is  balanced,  Vq  =  0,  the 
equation  reduces  to  become  ,        ,      ,  • 

-  (2.3) 


CRef        1  -  R 

If  Cfle/  is  known,  the  value  of  Cp  can  be  determined  using  the  balance  point  R.  However, 
the  value  of  Cp  is  not  needed  in  calculating  the  pressure  change  of  the  sample  because, 
through  a  calibration,  the  pressure  can  be  expressed  as  a  function  of  R,  as  discussed  in  a 
later  part  of  this  section. 

In  a  more  realistic  model,  it  can  be  shown  that  a  large  capacitance  in  the  output  cable 
can  significantly  reduce  the  bridge  sensitivity  [36].  This  problem  can  be  minimized  by 
using  a  short  length  of  cable.  Further  analysis  of  the  model  reveals  that  the  change  in 


'^Linseis,  Inc.,  P.  O.  Box  666,  Princeton  Jet.,  NJ  08550-0660,  phone:  (800)  732-6733. 
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out-of-phase  component  can  be  minimized  by  having  Cp  equal  to  CrsJ  and  similar  cable 
capacitance  values. 

To  minimize  heating  of  the  capacitor  plates  of  the  strain  gauge,  the  output  of  the 
function  generator  (HP3325A)  was  turned  on  only  when  a  measurement  was  made.  The 
HP3325A  was  used  instead  of  the  internal  oscillator  of  the  PAR124A  lock-in  amplifier 
because  its  output  could  be  controlled  by  a  computer.  The  on/off  timing  of  the  HP3325A 
was  remotely  controlled  by  a  LabWindow  program  which  was  run  on  the  316S  CompuAdd. 
The  program  was  written  by  Dr.  D.  A.  Sergatskov  and  was  designed  to  simultaneously 
control  three  HP  devices:  the  HP  3325A  Function  Generator,  the  HP  3478A  Multimeter, 
and  the  HP  3457A  Multimeter  (not  shown  in  Fig.  2.2). 

When  running,  the  program  first  sent  out  a  command  to  turn  on  the  output  of  the 
HP  3325A  function  generator.  The  output  signal  was  a  2. 12- Vpp  sine  wave  at  317  Hz  whose 
amplitude  and  frequency  could  be  set  in  the  program  menu.  The  frequency  was  chosen  to 
avoid  coupling  to  the  60  Hz  frequency  of  the  power  line.  With  an  excitation  of  2.12  Vpp, 
the  sensitivity  of  the  cell  was  approximately  5  Pa  (0.0007  psi)  at  11  MPa  (1600  psi).  Higher 
excitation  levels  increased  the  signal-to-noise  ratio  but  also  generated  more  heating  in  the 
plates,  which  was  not  desirable. 

Next  the  program  read  the  outputs  of  the  HP  3478A  and  HP  3457A  multimeters.  The 
HP  3478A  multimeter  was  used  exclusively  for  the  experimental  cell.  Depending  on  the 
stage  of  the  experiment,  the  HP  3457A  was  used  for  various  devices  such  as  the  melting 
pressure,  Ru02  resistance,  and  Pt  NMR  thermometers,  or  measurement  of  current  in  the 
sample  magnet.  A  typical  measurement  had  a  duration  of  90  seconds,  during  which  time 
the  program  read  and  stored  data  points  in  the  computer  memory  at  a  rate  of  2  points/sec. 
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Finally,  the  program  switched  oflF  the  HP  3325A  by  setting  the  output  voltage  level  to 
zero.  The  program  averaged  the  last  20  read  points  (10  seconds)  and  stored  the  value  in  a 
computer  file.  Only  these  last  points  were  reliable  because  of  the  30  second  time  constant, 
r,  set  on  the  PAR  124A.  The  lock-in  amplifier's  output  yields  correct  readings  only  after 
the  signal  has  been  turned  on  for  an  amount  of  time  greater  than  about  three  r. 

During  the  precooling  of  a  sample  (which  could  last  as  long  as  six  weeks),  one  mea- 
surement with  a  duration  of  90  seconds  was  made  every  hour  to  monitor  the  slow  change 
in  pressure.  In  a  34  hour  sample  demagnetization,  the  waiting  time  between  two  measure- 
ments was  reduced  from  1  hour  to  10  minutes,  and  then  to  5  minutes  for  fields  below  0.1  T 
in  order  to  obtain  better  resolution  of  the  pressure  change  of  the  sample  with  magnetic 
field.  Once  the  onset  of  ordering  was  observed,  data  were  taken  continuously  in  order  to 
provide  more  data  below  the  transition  point  since  the  magnetic  transition  occurs  at  very 
low  field.  With  a  demagnetization  rate  of  73  mT/hr,  the  sample  demagnetization  below 
the  transition  was  completed  in  a  few  minutes. 

2.3    Superconducting  Magnet 

In  this  work,  two  superconducting  magnets  were  made  with  maximum  fields  of  1  T 
and  2.5  T,  respectively,  in  vacuum. ^"^  The  main  difference  between  them  is  the  type  of 
superconducting  wire  used.  The  first  magnet  (1  T)  has  a  single  filament  superconducting 
wire,  and  the  second  (2.5  T)  a  multi-filament  wire,  the  construction  of  which  is  described 
in  the  following  text. 

'^The  1  T  magnet  was  initially  used  in  precooling  of  the  ^He  sample  but  was  replaced  by  the  2.5  T  magnet 
in  order  to  raise  the  precooling  temperature. 
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The  magnet  wire,  made  by  Supercon,  Inc./'*  has  54  superconducting  filaments  of 
52%  niobium-48%  titanium  embedded  in  a  0.257  mm  (0.0101")  dia.  copper  wire.  The  wire  is 
electrically  insulated  by  a  layer  of  Formvar  and  has  a  total  diameter  of  0.287  mm  (0.0113"). 

The  magnet  consists  of  two  parts  as  shown  in  Fig.  2.3  and  is  based  on  the  design 
of  Isareleson  and  Gould  [37].  The  first  part  consists  of  a  main  coil,  two  notches,  and  two 
vertical  compensation  coils.  The  notches  increase  the  uniformity  of  the  field  near  the  center 
of  the  magnet  and  the  vertical  compensation  coils,  wound  in  the  opposite  direction  of  the 
main  coil,  reduce  the  stray  field  along  the  z-axis.  The  suppression  of  the  vertical  field 
outside  the  magnet  is  particularly  important  because  an  aluminum  superconducting  heat 
switch  with  a  critical  field  of  only  11  mT  is  located  directly  above  the  magnet  and  any  large 
stray  field  would  reduce  its  effectiveness.  ' 

The  main  coil  consisted  of  21  layers,  each  having  264  turns.  The  notches  each  had  5 
layers  with  71  turns  per  layer  and  the  vertical  compensation  coils  had  14  layers  with  35 
turns  per  layer.  Upon  completion  each  layer  was  brushed  with  a  dilute  solution  of  GE 
varnish  and  acetone  to  prevent  it  from  loosening.  When  the  winding  of  the  first  part  was 
completed,  the  outer  layer  was  coated  with  Stycast  1266  epoxy.  The  whole  part  was  then 
immediately  placed  into  a  dessicator  and  pumped  to  remove  air  so  that  the  epoxy  could 
penetrate  into  the  inner  layers. 

The  second  part  is  an  active  shield  which  has  3  layers  (wound  in  the  opposite  direction 
of  the  main  coil)  with  372  turns  in  the  first  layer,  321  turns  in  the  second  layer,  and  270 
turns  in  the  third.  Like  the  first  part,  it  was  coated  with  dilute  GE  varnish,  and  epoxy  and 
an  additional  layer  of  dental  floss  was  added  to  provide  extra  protection  for  the  outermost 

'^Supercon,  Inc.,  830  Boston  Turnpike,  Shrewsbury,  MA  01545,  phone:  (508)  842-0174,  fax:  (508)  842- 
0847. 
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layer.  The  shield  was  made  in  an  attempt  to  cancel  out  the  stray  field  in  the  radial  direction, 
but  was  not  very  successful.  With  a  center  field  in  excess  of  0.73  T,  the  outside  field  was 
large  enough  to  prevent  the  nearby  Pt  NMR  thermometer  from  working.  Its  resonant  signal 
was  completely  destroyed  by  the  stray  field,  despite  having  a  superconducting  lead  shield 
surrounding  it.  With  a  0.73  T  center  field,  the  stray  field  at  the  thermometer  was  estimated 
to  be  20  mT.  However,  this  was  large  enough  to  interfere  with  the  14  mT  operating  field  of 
the  thermometer.  With  space  limitations  and  a  large  field  requirement  for  the  precooling 
of  the  sample,  the  problem  was  difficult  to  solve  and  meant  that  the  Pt  NMR  thermometer 
could  not  be  used  in  this  work.  However,  the  lack  of  the  Pt  NMR  thermometer  did  not 
create  a  significant  set  ba<:k  to  the  experiment  since  the  temperature  of  the  bundle  could  be 
measured  by  a  melting  pressure  thermometer  down  to  0.4  mK,'^  described  in  Section  2.5.1. 
Lower  temperatures  could  be  estimated  by  using  the  magnetic  field  on  the  Cu  bundle  and 
comparing  it  to  the  values  in  other  demagnetizations  in  which  the  Pt  NMR  thermometer 
had  been  used.  '    ••  .       f  '..:'),  . 

The  bore  of  the  magnet  has  a  diameter  of  2.84  cm  (1.12").  Its  large  diameter  is  de- 
signed not  only  to  work  with  the  current  pressure  cell  but  also  with  future  susceptibility 
experiments  which  will  require  a  large  space.  Including  the  a<;tive  shield,  the  total  size  of 
the  magnet  is  a  cylinder  of  diameter  8.1  cm  (3.2")  and  length  11.4  cm  (4.5").  The  magnet 
was  thermally  anchored  to  the  mixing  chamber  and  supported  by  three  Vespel-22  legs^^  to 
prevent  a  heat  leak  to  the  copper  stage.  The  legs  have  a  diameter  of  0.64  cm  (0.25")  and  a 
length  of  1.8  cm  (0.7").  Figure  2.4  shows  the  schematic  of  the  experimental  arrangement. 

'*The  lowest  temperature  needed  in  this  experiment  was  0.5  mK  for  100%  entropy  removal  in  a  sample. 
'^E.  I.  Du  Pont  De  Nemours  &  Company  (Inc.),  Du  Pont  Polymers,  Wilmington,  DE  19898,  phone:  (800)  - 
222-8377. 
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To  help  design  the  magnet  a  computer  program  written  in  FORTRAN  by  Chris  Gould 
at  the  University  of  Southern  California  was  used  to  calculate  the  field  profiles.  It  allows 
the  user  to  specify  the  dimensions  of  the  coils  and  their  locations  with  input  parameters 
of  the  coil  inner  radii,  outer  radii,  half  length  values,  vertical  offsets  to  the  main  coil,  and 
the  current  density.  Dennis  Greywall,  at  AT&T  Bell  Lab,  also  wrote  a  similar  program 
and  distributed  it  as  an  executable  file.  Except  for  its  inability  to  handle  the  vertical 
compensation  coils,  it  can  calculate  the  field  profile  for  a  magnet  which  consists  of  a  main 
coil,  notches,  and  an  active  shield.  The  Greywall  program  is  easier  to  use  and  was  employed 
in  designing  simple  magnets  such  as  that  of  the  heat  switch.  Figure  2.5  shows  how  well  the 
calculated  curve,  which  has  been  generated  by  Gould's  program,  agrees  with  data  along 
the  z-axis  of  the  magnet  measured  at  4  K  using  a  Hall  probe.  '• 

The  operation  of  the  superconducting  magnet  requires  a  persistent  switch  [38]  whose 
function  is  to  keep  a  current  in  the  superconducting  magnet  without  having  it  connected  to 
a  power  supply  at  all  times.  The  switch,  located  on  top  of  the  vacuum  can,  reduces  heating 
and  noise  dumped  into  the  cryostat.  The  construction  of  the  switch  is  quite  simple  and 
consists  of  three  main  parts:  a  stainless  steel  tube,  a  segment  of  superconducting  wire,  and 
a  heater  wire.  The  stainless  steel  tube  has  a  length  of  5  cm  (2")  with  a  0.64  cm  (0.25") 
diameter  and  a  wall  thickness  of  0.5  mm  (0.02").  The  superconducting  wire  is  multi-filament 
52%  niobium-48%  titanium,  the  same  type  of  wire  that  is  used  for  the  magnet.  Its  length  is 
roughly  1.2  m  (4  ft)  allowing  both  of  its  ends  to  reach  the  mixing  chamber  where  they  can 
be  clamped  together  with  the  wires  of  the  2.5  T  magnet.  The  Cu-Ni  heater  wire  is  61  cm 
long  and  0.08  mm  (0.003")  in  diameter  with  a  resistance  of  2.86  J2/cm  (87.08  fi/ft).  The 
heater  wire  was  tightly  wound  non-inductively  around  the  superconducting  wire  as  shown 
in  Fig.  2.6. 
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Figure  2.5.  Field  profile  of  the  2.5  Tesla  magnet. 
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The  superconducting  wire  was  bent  carefully,  to  prevent  making  a  kink,  at  its  mid-point. 
Both  wires  were  inserted  inside  the  stainless  steel  tube  and  Stycast  2850FT  epoxy  was  used 
to  fill  the  empty  space.  The  epoxy  keeps  the  wires  from  vibrating  and  provides  thermal 
insulation.  The  persistent  switch  requires  a  current  of  approximately  20  mA  to  keep  the 
superconducting  wire  in  the  normal  state. 

The  switch  was  tied  down  on  the  top  of  the  vacuum  can  with  the  superconducting  wire 
run  inside  the  can  down  to  the  mixing  chamber.  The  wire  was  carefully  thermal  anchored 
at  various  parts  of  the  refrigerator  (such  as  the  1  K  pot  and  still)  using  GE  varnish,  silver 
tape,  and  dental  floss.  At  the  mixing  chamber,  each  end  of  the  wire  was  connected  to  the 
magnet  leads  by  clamping  with  the  following  method. 

Before  the  wires  were  clamped,  about  5  cm  (2")  of  Formvar  insulation  was  stripped 
from  each  end  by  Strip-X.^^  The  copper  layer  was  then  etched  away  with  a  dilute  nitric 
acid  solution  until  the  superconducting  filaments  were  completely  exposed.  The  filaments 
from  the  magnet  and  the  persistent  switch  were  twisted  together  and  inserted  into  a  Cu  tube 
approximately  6  cm  (2.5")  in  length  and  0.5  cm  (0.2")  in  diameter  which  was  then  flattened 
using  a  press  with  a  pressure  of  660  MPa.  With  such  a  great  pressure,  an  impression  of  the 
filaments  could  be  seen  on  the  outside  of  the  flattened  copper  tubes.  It  was  important  to 
keep  the  pressure  below  1  GPa  to  prevent  stretching  of  the  copper  and  causing  the  filaments 
to  break.  Using  this  type  of  clamp  a  test  on  two  straight  wires  showed  that  the  filaments 
remained  superconducting  with  current  up  to  100  A  in  liquid  ''He.  The  clamping  of  the 
magnet  leads  worked  extremely  well:  after  six  weeks  of  persisting  at  2.5  T  (34  A)  during  a 
long  sample  precooling,  the  decay  of  current  in  the  magnet  was  not  detectable. 
'^GC  Electronics,  1801  Morgan  Street,  Rockford,  IL  61102. 
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To  energize  the  magnet,  the  heater  of  the  persistent  switch  is  turned  on  (step  A  of 
Fig.  2.7)  causing  the  superconducting  segment  to  become  normal.  The  normal  segment 
therefore  has  a  large  resistance  compared  with  the  superconducting  magnet  wire  so  that  all 
current  from  the  power  supply  flows  directly  to  the  magnet.  The  current  in  the  magnet  can 
then  be  increased,  and  when  the  maximum  current  is  reached,  the  heater  is  turned  off.  The 
wire  segment  becomes  superconducting  and  the  current  begins  to  flow  through  the  switch 
(step  B).  Initially,  the  currents  from  the  power  supply  and  the  magnet  are  equal,  so  the  net 
current  in  the  switch  is  zero.  However,  as  the  current  supply  is  ramped  down,  the  current 
that  flows  through  the  switch  increases  (step  C).  Finally  the  magnet  current  is  persisted 
when  there  is  no  current  remaining  in  the  power  supply  (step  D). 

To  remove  the  current  from  the  magnet,  the  steps  in  Fig.  2.7  are  carried  out  in  reverse 
order,  from  D  to  A.  The  current  inside  the  magnet  must  first  be  matched  by  that  of  the 
outside  power  supply.  The  current  inside  the  switch  gets  smaller,  step  C,  and  eventually 
becomes  zero  when  the  currents  have  the  same  value,  step  B.  The  heater  is  then  turned 
on,  step  A,  and  the  current  in  the  magnet  can  be  removed  by  reducing  the  current  in  the 
power  supply. 

The  magnet  was  energized  by  a  programable  100  A  Kepco  ATE6-100A  power  supply^® 
whose  output  current  can  be  controlled  remotely  by  applying  a  dc  voltage  between  0  and 
1  V  to  its  control  terminals.  A  ramp  unit  built  by  our  electronics  shop  automatically 
controls  the  ramping  process.  It  is  a  digital  linear  waveform  generator  which  consists  of 
two  main  parts.  The  first  is  a  microcontroller  board  which  allows  the  user  to  input  the 
target  current  and  ramp  time  via  a  computer  terminal  and  an  RS232  cable.  The  second 
part  is  a  16  bit  digital-to-analog  converter  (DAC)  board  whose  function  is  to  convert  the 

''Kepco,  Inc.,  131-38  Sanford  Avenue,  Flushing,  NY  11352,  phone:  (718)  461-7000. 
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Figure  2.7.  Steps  in  ramping  up  and  persisting  current  in  a  superconducting  magnet,  see 
text  for  more  detail. 
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digital  signal  from  the  microcontroller  to  an  analog  control  voltage  for  the  power  supply. 
The  full  output  range  of  the  ramp  unit  is  from  0  to  5  V  and  a  voltage  divider  is  used  to 
match  the  control  voltage  to  that  of  the  power  supply. 

2.4  Cryostat 

The  cryostat  is  suspended  from  an  aluminum  hexagonal  plate  with  side  length  0.70  m 
(27.75")  and  thickness  7.6  cm  (3")  consisting  of  three  layers.  The  top  and  bottom  layers  are 
made  of  solid  aluminum  and  each  has  a  thickness  of  1.3  cm  (0.5")  while  the  middle  layer, 
which  is  bonded  with  epoxy  to  the  outer  layers,  has  an  aluminum  honeycomb  structure  with 
thickness  5.1  cm  (2.0").  With  this  type  of  construction  the  plate  is  light  but  strong  enough 
to  support  the  weight  of  the  cryostat  without  bending.  The  plate  is  bolted  on  three  hollow 
aluminum  square  beams  which  are  welded  together  to  form  an  equilateral  triangle.  Three 
Vibraplane  1206  active-air  isolation  mounts,  made  by  Kinetic  Systems,  Inc.,^^  support  the 
triangle  at  its  corners.  The  mounts  are  set  on  top  of  a  4-m-(13')  tall  concrete  tripod  which 
sits  on  three  2-ton  concrete  blocks.  The  blocks  are  isolated  from  building  vibrations  by 
soft  rubber  gaskets.  This  arrangement  helps  to  minimize  mechanical  vibration  which  might 
couple  to  the  cryostat.  For  more  information,  see  Ref.  [39]. 

The  cryostat  can  be  divided  roughly  into  two  main  cooling  stages.  The  first  part  is 
the  dilution  refrigerator  and  the  second  part  is  the  copper  demagnetization  stage.  An 
aluminum  superconducting  heat  switch,  situated  between  the  two  parts,  controls  the  heat 
flow  between  them.  Additional  information  on  the  cryostat  can  also  be  found  in  Ref.  [40]. 

'^Kinetic  Systems,  Inc.,  20  Arboretum  Road,  P.O.  Box  K,  Roslindale,  MA  02131,  phone:  (617)  522-8700, 
fax:  (617)  522-6323. 
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2.4.1    Dilution  Refrigerator 

The  dilution  unit  is  an  Oxford  Instruments  model  1000.^°  The  volume  of  the  ^He/^He 
mixture  used  in  this  unit  is  approximately  250  NTP  liters  at  room  temperature  with  16.4% 
of '^He.-^^  The  mixture  is  stored  in  three  banks  of  storage  tanks  each  of  which  has  a  capacity 
of  117  1. 

The  gas  circulation  system  of  the  unit  is  composed  of  two  Edwards  9B3  diffusion 
pumps, ^•^  an  Edwards  EH250  mechanical  booster  (Roots),  and  an  Alcatel  2063H  rotary 
pump^"^  connected  in  series  in  the  cited  order.  The  9B3's  are  separated  from  the  still  line 
by  an  air-operated  gate  valve  and  a  freon-cooled  baffle  which  prevents  mist  of  the  pump  oil 
from  getting  into  the  cryostat. 

Figure  2.8  shows  the  gas  handling  system  of  the  cryostat  where  dashed  lines  represent 
hidden  gas  lines.  Valve  #10  is  a  relief  solenoid  valve  which  automatically  opens  when  the 
line  pressure  exceeds  90  KPa  (650  Torr).  Valve  #32  is  a  one-way  safety  valve  and  #23  and 
#36  are  micro-valves.  When  the  dilution  unit  is  running,  the  ^He  with  about  10%  ''He  is 
pumped  from  the  still  (F)  by  the  9B3's.  The  mixture  circulates  through  the  Roots  and 
Alcatel  pumps  and  returns  to  the  gas  panel  from  the  back  side  of  the  Alcatel.  Before  going 
back  to  the  cryostat,  the  mixture  goes  through  one  of  the  two  LN2  cold  traps,  a  Hasting 
U5KM  flow  meter.^-*  and  valve  #29. 

^"Oxford  Instruments,  130A  Baker  Avenue,  Concord,  MA  01742,  phone:  (508)  369-9933,  fax:  (508)  369- 
6616. 

^'This  data  is  obtained  from  measurement  made  on  3/2/1990  during  testing  #3  of  the  cryostat. 

Edwards  High  Vacuum  International,  One  Edwards  Park,  301  Ballardvale  Street,  Wilmington,  MA 
01887,  phone:  (800)  848-9800,  fax:  (508)  658-7969. 

Alcatel,  40  Pond  Park  Road,  South  Shore  Park,  Hingham,  MA  02043,  phone:  (617)  749-8710,  fax:  (617)- 
749-8660. 

^^Teledyne  Brown  Engineering,  Hasting  Instruments,  P.  O.  Box  1436,  Hampton,  VA  23661, 
phone:  (800)  950-2468,  fax:  (804)  723-3925. 
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Figure  2.9.  Cooling  power  of  the  Oxford  1000  dilution  unit. 

Figure  2.9  shows  the  cooling  power  of  the  unit  from  110  mK  down  to  5  mK  taken  from 
measurements  made  in  October  of  1992  after  the  removal  of  some  coaxial  cables  which  had 
caused  a  large  heat  leak  to  the  mixing  chamber.  The  circulation  rate  is  that  of  the  whole 
mixture,  measuring  ^He  and  '*He  combined.  The  cooling  power  of  the  mixing  chamber  can 
be  expressed  by  [41] 

dQIdt  =  aTl,  -  b  (2.4) 

where  Tmc  is  the  temperature  of  the  mixing  chamber  and  a  and  6  the  fitting  parameters, 
except  near  the  minimum  temperature  where  the  fit  is  not  valid. 
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2.4.2    Aluminum  Superconducting  Heat  Switch 

The  copper  bundle  can  be  thermally  isolated  from  the  dilution  unit  with  an  aluminum 
superconducting  heat  switch.  Aluminum  becomes  superconducting  at  about  1.2  K  in  zero 
field  and  has  a  critical  field  Bc(0)  of  about  11  mT.  The  thermal  conductivity  of  aluminum 
is  electronic  and  is  thus  linear  with  temperature  in  the  normal  state  (k„),  with  a  de- 
pendence when  it  becomes  superconducting  (kj).  For  T  <  1.2  K,  the  switch  can  be  closed 
(normal)  by  applying  a  field  greater  than  Be  or  opened  (superconducting)  by  removing  the 
applied  field.  In  1978,  Mueller  et  al.  [42]  constructed  the  first  aluminum  heat  switch.  Their 
switch  has  a  large  switching  ratio  k„/k5  of  1600/T'^  at  58  mK  compared  to  the  typical  val- 
ues of  lOO/T^  to  300/T^  of  a  more  popular  tin  (Sn)  switch.  One  main  problem  encountered 
by  Mueller  et  al.  in  making  an  aluminum  heat  switch  was  that  its  surface  becomes  oxidized 
very  quickly  when  exposed  to  air.  This  problem  was  solved  by  electroplating  the  aluminum 
with  gold,  a  procedure  described  in  Ref.  [42]. 

In  precooling  of  the  copper  bundle  the  cooling  rate  depends  on  the  thermal  conductance 
between  the  mixing  chamber  and  the  bundle.  For  this  cryostat  we  found  that  conductance 
between  the  contact  of  the  heat  switch  and  the  thread  annulus  of  the  mixing  chamber  was 
poor  (105r  mW/K^),  a  serious  problem  in  the  design  of  the  mixing  chamber  by  Oxford. 
This  problem  was  overcome  by  constructing  a  heat  switch  which  can  be  bolted  to  24  avail- 
able tapped  holes  at  the  bottom  of  the  mixing  chamber.  This  arrangement  provided  a 
measured  conductance  of  540r  mW/K^,  an  improvement  over  Oxford's  design  by  a  factor 
of  5. 

The  conductance  of  the  heat  switch  itself  was  also  improved  from  210  to  390r  mW/K^ 
by  increasing  the  thickness  of  the  aluminum  strips  from  0.5  mm  to  1  mm.  With  these  two 
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improvements,  the  net  conductance  between  the  mixing  chamber  and  the  copper  bundle 
was  increased  from  61  to  200T  mW/K^  which  reduces  the  cooling  time  by  half  [40]. 

2.4.3    Copper  Demagnetization  Stage 

The  stage  is  made  from  350  moles  of  Cu  consisting  of  19  Cu  plates,^^  each  of  which 
has  a  thickness  of  3.2  mm  (0.125").  The  plates  have  a  length  of  91  cm  (36")  and  a  width 
varying  from  11  mm  (0.43")  to  63  mm  (2.5").  When  put  together,  the  plates  have  the 
form  of  a  cylinder  of  63  mm  (2.5")  in  diameter,  see  Fig.  2.10.  The  plates  are  electron-beam 
welded  to  copper  flanges  at  each  end.  The  flanges  are  gold  plated  to  reduce  oxidation  and 
provide  good  thermal  contact  for  attaching  experiments. 

The  performance  of  the  dilution  unit  during  a  precooling  of  the  copper  nuclear  demag- 
netization stage  is  shown  in  Fig.  2.11.  The  expression  for  the  precooling  time  t  required  for 
the  stage  to  reach  a  certain  temperature  T  can  be  derived  in  the  following  way.  First,  the 
heat  flow  from  the  copper  stage  to  the  mixing  chamber  is  given  by 


dQ/dt='^{T'-TL).      '  '  '  (2.5) 


Using  Eq.  2.4  to  eliminate  T^^,  one  obtains 


k 

dQ/dt=  j^{aT^  -b).  (2.6) 
2  -|-  a 


The  nuclear  heat  capacity  of  the  copper  stage  has  the  form 


C=Y2  (2-7) 


^*With  173  of  the  350  moles  is  in  a  8  T  magnet.  The  use  of  many  plates  in  place  of  a  single  piece  is 
to  minimize  the  cross-section  area  of  the  copper  stage  and  thus  reduce  the  eddy-current  heating  during  a 
demagnetization. 
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Figure  2.10.  Cross-section  view  of  the  copper  bundle. 
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Figure  2.11.  Temperature  versus  time  during  the  precooling  of  the  copper  demagnetization 
stage. 
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which  implies  that  the  cooling  (minus  sign)  of  the  stage  can  be  expressed  as 


dQ/dt=-^dTldt.  (2.8) 


Now,  substituting  this  into  Eq.  2.6  and  integrating,  we  have 
Performing  the  integration  yields  [40] 


where  To  =  \/o/a  is  the  minimum  temperature  obtained  by  setting  dQ/dt  to  zero  in  Eq.  2.4. 
The  values  for  the  theoretical  line  in  Fig.  2.11  are  k  =  130  mW/K^,  a  =  SlmW/K^, 
b  —  1.7  //W,  and  c  =  25  mJK  all  of  which  were  obtained  from  direct  measurements  or  a 
calculation.  In  a  field  of  7.02  T  it  takes  about  72  hours  to  precooling  the  copper  bundle  to 
11  mK. 

The  demagnetization  schedule  of  the  copper  stage  has  the  following  steps: 
7.02  T        3.74  T       1.89  T  ^-H'  0.93  T  ^•^'0.46  T  H^'  0.23  T         0.19  T       0.93  T 
"^■H-  0.47  T  "^-H-  0.23  T  '^■H'  0.12  T  H^'  0.0058  T. 

This  results  in  the  demagnetization  profile  of  the  bundle  shown  in  Fig.  2.12.  The  straight 
line  is  for  an  ideal  adiabatic  demagnetization  given  by 


T  =  r,^,  (2.11) 
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Figure  2.12.  Demagnetization  curve  of  the  copper  stage. 
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where  T,  and  B,-  are  the  initial  temperature  and  magnetic  field  respectively,  and  B  is 
the  current  field.  The  measured  curve  roughly  follows  the  ideal  curve  to  approximately 
1  mK.  Below  this  temperature,  the  deviation  from  ideality  becomes  larger  due  to  parasitic 
heat  capacities  and  external  heat  leaks.  The  minimum  measured  temperature  achieved 
by  this  copper  stage  is  45  nK.^^  Figure  2.12  has  two  data  points  in  the  lowest  field.  The 
high  temperature  point  was  obtained  at  the  end  of  the  demagnetization  process  and  the 
low  point  after  closing  the  fill  line  of  the  1  K  pot.  Closing  the  line  helps  to  reduce  the 
vibrational  heating  and  allows  the  thermometer  to  reach  a  lower  temperature. 

2.5  Thermometry 

2.5.1    Melting  Pressure  Thermometer 

The  melting  pressure  thermometer  and  electronic  configuration  are  similar  to  those  of 
the  experimental  cell.  However,  because  it  operates  in  zero  field  and  at  a  lower  pressure, 
3.45  MPa  (500  psi),  the  physical  properties  of  the  thermometer  are  slightly  different  from 
those  of  the  experimental  cell.  The  strain  gauge  is  made  of  beryllium  copper  (BeCu)  with 
a  diaphragm  of  diameter  1.27  cm  (0.5")  and  thickness  0.74  mm  (0.029").  The  movable 
plate  is  epoxied  onto  a  post  at  the  center  of  the  diaphragm,  using  Stycast  2850FT.  The 
fixed  plate  is  also  fixed  to  the  holder  with  the  same  epoxy.  The  plates  have  a  diameter  of 
1.14  cm  (0.45")  and  are  separated  by  a  mylar  spacer  of  thickness  0.05  mm  (0.002").  The 
holders  are  held  together  by  four  2-56  stainless  steel  screws. 

The  body  of  the  thermometer  is  machined  out  of  a  99.999%  pure  silver  rod  with  a 
diameter  of  1.3  cm  (0.5")  and  a  length  of  1.9  cm  (0.75").  A  coin-silver  flange  is  soldered 
with  eutectic  silver-copper  solder  (71.9%  Ag  -  28.1%  Cu)  onto  the  rod  at  one  end.  The 

^^This  temperature  was  measured  by  a  Pt  NMR  thermometer  located  on  top  of  the  copper  bundle.  The 
temperature  at  the  center  of  the  bundle  is  assumed  to  be  lower  thcin  this  value. 
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flange  has  eight  tapped  holes  where  the  strain  gauge  is  bolted.  An  indium  0-ring  is  used 
to  provide  a  vacuum  seal  between  the  body  and  the  strain  gauge.  The  other  end  of  the 
rod  is  welded  to  a  silver  plate  which  has  eight  clearance  holes  for  screws  that  are  used  to 
attach  it  to  the  top  of  the  copper  demagnetization  stage.  Thermal  contact  between  the 
cell  body  and  the  -^He  is  provided  by  silver  (Ag)  sinter  and  eight  silver  rods  of  diameter 
0.8  mm  (0.03"). 

The  excitation  for  capacitance  measurements  is  provided  by  the  internal  oscillator  of 
an  Ithaco  393  lock-in  amplifier.  Because  the  thermometer  is  operated  at  melting  density, 
the  liquid  provides  an  excellent  thermal  path  to  the  silver  sinter;  so  the  excitation  voltage 
can  be  kept  on  at  all  times  without  concern  about  local  heating.  Also,  the  stainless  steel 
screws  used  on  the  thermometer  help  to  conduct  heat  away  from  the  ^He.  Note  that  in 
the  experimental  cell,  there  are  no  screws  nor  liquid  to  remove  the  heat  and  it  is  therefore 
necessary  to  turn  off  the  excitation  after  a  measurement  to  prevent  heating  of  the  sample 
in  the  experimental  cell. 

The  operating  frequency  for  the  thermometer  is  1.1  kHz  with  an  amplitude  of  1.22  Yrms- 
A  one-to-one  isolation  transformer  is  used  to  keep  the  system  from  shorting  to  ground.  The 
thermometer  and  its  reference  capacitor  (which  sits  on  the  mixing  chamber  of  the  dilution 
refrigerator)  are  connected  to  a  General  Radio  1493  ratio  transformer  via  three  coaxial 
cables  as  in  the  experimental  cell  configuration.  The  output  signal  from  the  bridge  is  fed 
into  the  lock-in  amplifier  detected,  amplified,  and  then  recorded  by  a  chart  recorder  and  a 
multimeter  (HP  3457A)  which  is  connected  to  a  personal  computer. 

Similar  to  that  of  the  experimental  cell,  the  strain  gauge  of  the  thermometer  is  trained 
and  calibrated  at  high  temperature.  See  Sec.  3.1.2  for  detailed  description  of  the  procedures. 
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The  temperature  can  be  calculated  from  the  pressure  using  the  fitting  parameters  given 
in  Ref.  [9].  For  temperature  higher  than  T/v  =  0.934  mK  (the  Neel  temperature  at  which 
the  ^He  solid  orders  into  the  LFP),  the  pressure  function  has  the  form 

5 

P{T)  -Pn=  Yl  ^n^"'  (2-12) 

n=— 4 

where  Pn  is  the  pressure  measured  at  T/v,  the  reference  point.  The  A„'s  are  the  fitting 

coefficients  and  have  the  following  values  for  P  in  kPa  and  T  in  mK: 

A_4  =  0.78405848  A_3  =  -4.5968629  A_2  =  9.5436359 

A_i  =  -10.041752  Ao  =  8.5580332  Ai  =  -4.4431076 

A2  =  1.5962639  x  lO'^        A3  =  -3.8815407  x  10"^    A4  =  7.1230765  x  10"^ 

A5  =  -6.0945626  x  IQ-^^. 

For  temperatures  below  Tyv,  the  temperature  can  be  determined  from  the  expression 

^      P{T)-Pn  =  Ao-AiT\       -  ,  (2.13) 

r    ••   ■  ■ 

with  Ao  =  0.1987  kPa  and  Ai  =  0.2611  kPa/mK*.        '  '  ;       ■  .  ;  ;  <  ' 

2.5.2    Resistance  Thermometers 

There  are  total  of  17  resistance  thermometers  at  different  locations  in  the  cryostat. 
Six  of  them  are  Speer  thermometers  provided  by  Oxford,  the  manufacturer  of  the  dilution 
unit.  The  others  (R7  -  R22)  are  RUO2  thermometers  [43]  which  have  been  calibrated  and 
provided  by  Professor  Mark  Meisel.  Except  for  R7  and  R8,  which  are  1  kfi  resistors,  the 
rest  are  500  f2.  The  resistance  values  of  these  thermometers  at  different  temperatures  are 


60 


listed  in  Table  2.1  to  Table  2.3  as  determined  by  using  a  Linear  Research  ac  resistance 
bridge  (LRllO)." 

In  this  experiment  the  most  important  thermometers  are  R4  and  R16  which  are  located 
on  the  mixing  chamber.  The  resistor  R4  is  used  in  forming  and  annealing  the  high  density 
^He  samples,  procedures  which  are  discussed  in  the  next  chapter.  The  resistor  R16  is 
used  to  monitor  the  temperature  of  the  dilution  refrigerator  after  the  heat  switch  has  been 
opened. 


Linear  Research,  Inc.,  5231  Cushman  Place,  Suite  21,  San  Diego,  CA  92110,  phone:(619)299-0719. 


61 


Table  2.1.  List  of  the  resistance  values  of  the  Oxford  refrigerator's  Speer  thermometers  at 
various  temperatures.  


Labels 

Rl 

R2 

R3 

R4 

R5 

R6 

Type  (Q) 

470 

470 

100 

470 

100 

220 

Locations 

Film 
Burner 

Still 

Cold 
Plate 

Mix. 
Chamber 

Mix. 
Chamber 

1  K 
Plate 

300  K  (kfi) 

0.80 

0.80 

0.43 

0.89 

0.46 

0.53 

77  K  (kQ) 

0.86 

0.87 

0.43 

0.98 

0.47 

0.55 

4.2  K  (kf2) 

1.20 

1.23 

0.46 

1.42 

0.49 

0.69 

Table  2.2.  List  of  the  resistance  values  of  the  RuOa  thermometers  above  the  heat  switch 
at  various  temperatures.  


Labels 

R7 

R8 

R14 

R15 

R16 

R17 

Locations 

Still 

Cold 

6  th  Heat 

7  th  Heat 

Mix. 

Heat  Switch 

Shield 

Shield 

Exchanger. 

Exchanger 

Chamber 

Top 

300  K  (kQ) 

1.31 

1.32 

0.97 

0.97 

0.94 

0.96 

77  K  (kfi) 

1.29 

1.30 

0.94 

0.94 

0.92 

0.94 

4.2  K  (kQ) 

1.61 

1.64 

0.89 

0.88 

0.89 

0.91 

4  mK(kf2) 

14.4 

19.4 

Table  2.3.  List  of  the  resistance  values  of  the  Ru02  thermometers  below  the  heat  switch 
at  various  temperatures.  


Labels 

R18 

R19 

R20 

R21 

R22 

Locations 

Heat  Switch 

Cage 

Cage 

Bundle 

Bundle 

Bottom 

Top 

Bottom 

Top 

Bottom 

300  K  (kfi) 

0.96 

0.98 

0.98 

0.98 

1.07 

77  K  (kQ) 

0.94 

0.96 

0.96 

0.96 

1.04 

4.2  K  (kQ) 

0.92 

0.93 

0.91 

0.91 

0.92 

4  mK(kQ) 

18.5 

18.8 

17.0 

27.7 

25.2 

CHAPTER  3 
EXPERIMENTAL  PROCEDURES 

In  this  chapter,  the  calibration  procedure  of  the  strain  gauges  of  the  experimental  cell 
and  melting  pressure  thermometer  is  covered  in  detail,  followed  by  the  discussion  on  the 
magnetic-field-,  ramp-rate-,  and  temperature-dependence  of  the  experimental  cell.  A  short 
section  on  temperature  regulation  method  of  the  cryostat  is  given  next,  proceeded  by  a 
detailed  description  of  forming  a  high-density  "^He  sample. 

3.1    Calibration  of  the  Strain  Gauges 

The  first  half  of  this  section  is  devoted  to  the  training  and  calibration  of  the  experimental 
cell  and  the  second  half  is  for  the  melting  pressure  thermometer. 

3.1.1    Experimental  Cell 

In  order  to  translate  the  capacitance  of  the  strain  gauge  to  the  pressure  of  the  sample, 
the  gauge  must  be  accurately  calibrated.  For  high  density  solid  ^He,  the  pressure  calibration 
was  carried  out  at  4  K,  above  the  solidification  point. 

The  first  step  of  the  calibration  process  is  to  train  the  diaphragm  of  the  gauge  to 
reduce  any  hysteresis  which  may  exist  in  the  diaphragm.  The  training  is  done  by  slowly 
pressurizing  the  cell  with  ^He  gas  from  0.0  to  14.5  MPa  (2100  psi)  in  15  minutes.  This 
process  is  carried  out  by  inserting  the  ''He-temperature  dipstick  of  the  ^He  gas  handling 
system,  shown  in  Fig.  3.1,  into  a  liquid  helium-4  (L^He)  dewar.  As  the  stick  cools,  it 
cryo-pumps  ^He  gas  from  the  storage  bottle  through  the  valves  Nl,  #1,  and  #2,  where 
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"N"  denotes  a  Nupro  valve^  and  all  other  valves  are  HIP  high  pressure  valves.^  The  thicker 
lines  represent  0.64  cm  (0.25")  o.d.  Cu  tubes  and  the  thinner  lines  1.5  mm  (0.06")  o.d. 
stainless  steel  capillaries.  The  gas  is  passed  through  a  LN2  cold  trap  to  remove  impurities 
such  as  air.  When  sufficient  gas  has  been  removed  from  the  bottle,  valve  #1  is  closed 
and  the  dipstick  slowly  pulled  out  of  the  L'*He  dewar  increasing  the  pressure  of  "^He  gas 
as  the  stick  warms  up.  With  valve  #3  opened,  the  gas  pressure  can  be  monitored  by  a 
2000  psi  (14  MPa)  Heise  gauge.  Before  entering  the  cryostat  through  valve  #4,  the  ^He 
gas  is  passed  through  another  cold  trap  to  ensure  that  there  is  no  air  contamination.  After 
reeiching  a  pressure  of  15  MPa  (2100  psi),  the  gas  is  slowly  removed  by  inserting  the  dipstick 
back  into  the  dewar.  When  the  pressure  drops  to  7  MPa  (1000  psi),  the  cell  is  then  reloaded 
with  ^He  gas  at  the  same  rate  to  14  MPa  (2000  psi).  From  our  experience,  this  was  the 
most  efficient  way  of  training  the  diaphragm. 

Step  two  of  the  process  is  the  calibration  itself  which  is  carried  out  by  maintaining 
the  cell  at  different  values  of  pressure  and  then  balancing  the  capacitance  bridge  to  deter- 
mine the  corresponding  capacitance  ratio.  The  calibration  is  done  against  the  Heise  gauge 
which  has  a  resolution  of  7  kPa  (1  psi).  The  calibration  range  for  this  cell  is  between 
8  MPa(1200  psi)  and  14  MPa  (2000  psi).  From  14  MPa  (2000  psi)  the  pressure  is  lowered 
in  340  kPa  (50  psi)  steps  until  it  reaches  11  MPa  (1600  psi).  Then  from  this  pressure  down 
to  10  MPa  (1400  psi),  the  step  size  is  reduced  to  70  kPa  (10  psi)  to  give  a  more  accurate 
calibration  in  the  low  temperature  working  range  of  the  cell.  Between  10  MPa  (1400  psi) 
and  8  MPa  (1200  psi),  the  pressure  step  is  increased  back  to  340  kPa  (50  psi).  When  the 
minimum  pressure  of  8  MPa  (1200  psi)  is  reached,  the  cell  is  reloaded  with  ^He  gas  so 

'Jax  Valve  k  Fitting  Company,  3633  Southside  Blvd.,  P.  O.  Box  16097,  Jacksonville,  FL  32245, 
phone:  (904)  642-2093. 

^High  Pressure  Equipment  Company,  1222  Linden  Avenue,  Erie,  PA  16505-3522,  phone:  (814)  838-2028, 
fax:  (814)-838-6075. 
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Figure  3.1.  ^He  gas  handling  system. 
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that  the  calibration  with  increasing  pressure  can  be  obtained.  The  pressure  steps  are  the 
same  as  in  the  calibration  for  decreasing  pressure.  Upon  reaching  the  maximum  pressure  of 
14  MPa  (2000  psi)  the  gas  is  again  removed  in  steps  for  another  calibration  with  decreasing 
pressure.  In  a  typical  calibration,  there  are  two  "down"  and  one  "up"  pressure  calibration 
processes.  By  plotting  the  pressure  versus  capacitance-ratio-reading  for  all  three  processes 
one  can  determine  the  hysteresis  of  the  cell. 

Figure  3.2  shows  the  calibration  curves  in  which  circles  and  triangles  denote  calibra- 
tions with  decreasing  pressure,  and  squares  the  calibration  with  increasing  pressure.  A 
polynomial  equation  of  the  form 

P  =  ao  +  aiR  +  a2R^ +  a3R^,  (3.1) 

was  used  to  fit  the  down-pressure  data  points,  where  P  is  the  value  of  the  pressure,  R 
the  corresponding  ratio  transformer  reading,  and  the  a,'s  are  fitting  parameters.  The 
hysteresis  of  the  cell  can  be  determined  by  subtracting  the  fitting  polynomial  curve  from 
the  up  pressure  data,  and  this  difference  is  shown  in  Fig.  3.3.  Although  this  cell  has  a 
hysteresis  of  about  30  kPa  (4  psi)  for  the  calibration  range  between  8  MPa  and  14  MPa,^ 
in  practice  the  hysteresis  is  much  smaller  because  of  a  narrower  working  pressure  range  of 
the  cell  at  low  temperature.  From  the  temperature  at  which  the  sample  becomes  solid,  the 
total  pressure  change  down  to  absolute  zero  is  approximately  0.34  MPa  (50  psi),  16  times 
less  than  the  calibration  range.  Furthermore,  during  the  precooling  and  demagnetization 
of  the  sample  the  total  pressure  change  is  less  than  400  Pa  making  the  hysteresis  during  a 
measurement  much  smaller  than  the  30  kPa  measured  during  the  calibration. 
^It  is  possible  that  the  Heise  pressure  gauge  itself  has  hysteresis. 
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Figure  3.2.  Calibration  curve  of  the  experimental  cell. 
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Figure  3.3.  Residuals  of  the  calibration  curve  of  the  experimental  cell. 
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3.1.2    Melting  Pressure  Thermometer 

Prior  to  the  calibration  at  around  1  K,  the  diaphragm  is  trained.  During  the  process, 
■'He  is  introduced  into  the  thermometer  from  zero  pressure  to  3.7  MPa  (530  psi)  in  5  minutes 
and  then  removed  at  the  same  rate.  When  all  the  gas  is  out,  it  is  loaded  back  into  the 
system  at  a  slower  rate,  about  3.4  kPa/sec  (0.5  psi/sec).  The  pressure  is  monitored  by  two 
pressure  gauges,  a  Heise  gauge  and  a  Paroscientific  pressure  transducer.^  The  transducer 
has  a  sensitivity  of  3.4  Pa  (0.0005  psi)  but  can  only  withstand  a  maximum  pressure  of 
6.9  MPa  (1000  psi)  (which  makes  it  unsuitable  for  calibration  of  the  experimental  cell.) 
When  the  pressure  reaches  3.6  MPa  (520  psi),  the  training  is  completed. 

The  calibration  range  is  between  3.6  MPa  (520  psi)  and  2.8  MPa  (400  psi)  and  is  initiated 
by  removing  the  gas  in  small  pressure  steps  of  70  kPa  (10  psi)  down  to  3.5  MPa  (500  psi). 
At  every  step,  readings  are  taken  only  when  the  pressure  change  in  2  seconds  is  less  than 
690  Pa  (0.01  psi),  as  determined  from  the  Paroscientific  gauge.  Between  3.5  MPa  (500  psi) 
and  3.4  MPa  (490  psi),  the  pressure  steps  are  reduced  to  7  kPa  (1  psi)  because  this  is  the 
operating  range  of  the  thermometer  at  low  temperature.  Below  3.4  MPa  (490  psi),  the 
steps  are  increased  back  to  70  kPa  (10  psi)  until  2.8  MPa  (400  psi).  Like  the  calibration 
of  the  experimental  cell,  there  are  two  calibrations  with  decreasing  pressure  and  one  with 
increasing  pressure. 

In  the  next  section,  we  discuss  the  magnetic-field-,  ramp-rate-,  and  temperature-dependence 
of  the  experimental  cell.  These  effects  must  be  subtracted  from  the  raw  pressure  data  to 
obtain  the  true  pressure  change  in  the  sample  during  the  precooling  and  demagnetization. 
^Paroscientific,  4500  NE  148  th  Avenue  Redmond,  WA  98052-5126,  phone:  (206)  883-8700. 
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3.2    Magnetic-Field-Dependence  of  the  Cell 

In  our  first  experimental  cell,  we  encountered  a  large  magnetic-field-dependence  which 
was  mainly  created  by  four  stainless  steel  screws  used  to  hold  the  strain  gauge  together. 
Figure  3.4  shows  the  field-dependence  of  the  first  cell  expressed  as  an  equivalent  pressure 
change  of  more  than  1.6  kPa  between  0  and  1.6  T.  This  huge  field-dependence  is  about  4 
times  larger  than  the  expected  400  Pa  change  in  the  bcc  solid  pressure  during  a  demag- 
netization and  made  observation  of  the  actual  pressure  change  in  the  sample  difficult.  For 
this  reason  a  second  cell  was  made  in  an  attempt  to  reduce  the  effect  of  the  magnetic  field. 

The  second  sample  cell  has  no  stainless  steel  screws  and  exhibits  a  much  smaller  field- 
dependence,  as  shown  in  Fig.  3.4.  The  field-dependence  is  larger  at  high  fields  with  a 
total  change  of  270  Pa  between  1.8  and  2.5  T  and  the  equivalent  pressure  increases  with 
decreasing  field.  For  fields  from  1.8  T  down  to  0  T,  the  total  change  is  less  than  70  Pa 
(approximately  4%  of  that  in  the  first  cell)  and  reverse  in  direction  as  shown  in  Fig.  3.5. 
Fortunately,  the  dependence  in  low  fields  was  reproducible  and  the  pressure  change  in  the 
sample  can  be  obtained  by  subtracting  the  field-dependence  from  the  raw  pressure  data. 
The  field-dependence  was  fitted  to  a  polynomial  form  which  was  used  in  analyzing  all 
sample  demagnetization  data.  The  fitting  form  for  low  fields  (0  <      <  1.8  T)  is 

AP  =  -2.41925  +  67.18085o  -  7.158285^  -  4.35145^,  (3.2) 

with  AF  in  units  of  Pascals  and  Bo  in  Tesla.  The  high-field  portion  (1.8  <  So  <  2.5  T)  can 
also  be  fitted  with  a  polynomial  form,  though  this  portion  varies  from  one  experimental 
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B  (T) 

Figure  3.4.  Magnetic-field-dependence  in  decreasing  field  of  the  first  and  second  experimen- 
tal cells  at  fixed  temperature.  The  first  cell  has  four  stainless  steel  screws  and  the  second 
cell  has  none. 
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Figure  3.5.  Magnetic-field-dependence  of  tiie  second  experimental  cell  at  30  mK  in  decreas- 
ing fields  .  The  empty-cell  curve  was  measured  when  there  was  no  sample  in  the  cell  while 
the  filled-curve  was  obtained  with  a  solid  sample.  A  slight  diflference  in  the  two  curves  may 
be  due  to  the  diflferent  shapes  of  the  diaphragm  with  and  without  sample. 

run  to  another.  For  the  filled-cell  data  in  Fig.  3.5,  it  is  given  by 

AF  =  2144.74  -  3678.045o  +  2188.865^  -  436.904B^.  (3.3) 

The  cell  also  exhibits  a  field  ramp-rate-dependence  as  shown  in  Fig.  3.6.  A  jump 
in  the  equivalent  pressure  occurred  whenever  there  was  a  change  in  the  ramp-rate.  The 
equivalent  pressure  increases  when  the  ramp-rate  was  decreased  and  decreases  when  the 
rate  was  increased.  The  size  of  the  jumps  depended  on  the  ramp-rate  and  the  field  at  which 
it  occurs.  The  size  was  smaller  in  high  fields  or  when  the  change  in  the  ramp-rates  was 
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Figure  3.6.  Ramp-rate-dependence  of  the  experimental  cell  at  fixed  temperatures. 
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small.  The  upper  curve  in  Fig  3.6  was  obtained  at  150  fiK  with  an  initial  ramp-rate  of 
-73  mT/hr.  At  the  field  of  approximately  0.17  T,  the  ramp-rate  was  reduced  to  -6  mT/hr 
which  caused  a  jump  of  approximately  60  Pa.  When  the  ramping  stopped  at  zero  field, 
a  jump  of  about  80  Pa  was  measured  during  5  hours  of  waiting.  The  bottom  curve  was 
obtained  at  4  mK  with  an  initial  ramp-rate  of  -146  mT/hr.  At  the  fields  of  1.38  and  0.69  T, 
the  ramping  was  stopped  for  1  hour  and  two  equivalent  pressure  increases  of  40  and  50  Pa, 
respectively,  were  observed.  At  the  field  of  0.14  T,  the  ramp-rate  was  reduced  to  -37  mT/hr 
which  caused  a  jump  of  approximately  20  Pa.  When  the  ramping  was  stopped  at  zero  field, 
the  equivalent  pressure  increased  by  70  Pa  over  4  hours  of  waiting.  The  measurements 
of  the  ramp-dependence  were  made  in  an  empty  cell,  but  this  effect  could  still  be  seen  in 
the  filled  cell.  In  order  to  avoid  this  dependence  during  a  sample  demagnetization,  the 
demagnetization  rate  was  kept  constant  at  all  times.  Since  there  existed  a  large  equivalent 
pressure  jump  at  zero  field,  all  data  taken  after  the  demagnetization  were  not  included  in 
the  data  analysis.        '  •  ' 

We  do  not  know  the  cause  of  this  ramping  effect  but  believe  that  it  is  not  due  to 
eddy  current  heating.^  In  Fig.  3.7  we  have  plotted  the  square  of  the  ramp-rate  versus  the 
difference  in  the  squares  of  the  temperatures  of  the  capacitance  plates  and  those  of  the 
copper  bundle.  The  temperature  of  the  plates  was  determined  by  calculating  the  size  of  the 
pressure  jumps  and  relating  these  values  to  the  temperature-dependence  (discussed  in  the 
next  section)  of  the  cell.  If  the  ramping  effect  were  caused  by  the  eddy-current  heating,  the 
fitted  curve  would  have      proportional  to  T^^^^^  —  T^^^^j^,  but  we  find  no  such  dependence. 


An  induced  magnetic  force  on  the  diaphragm  can  also  be  discounted  because  this  force  would  cause  an 
increase  in  the  equivalent  pressure  with  a  magnitude  of  only  10~^  to  what  observed. 


Figure  3.7.  Ramp-rate  versus  temperature.  The  fit  line  shows  that  the  heating,  if  real,  is 
not  from  the  metal  plates. 
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Figure  3.8.  Temperature-dependence  of  the  experimental  cell  when  it  is  empty. 

3.3  Temperature-Dependence 

Besides  the  field  and  ramp-rate-dependence,  the  cell  also  has  a  small  temperature- 
dependence.  The  temperature-dependence  was  measured  for  an  empty  cell^  and  is  shown 
in  Fig.  3.8.  For  temperature  below  10  mK,  data  were  taken  in  only  zero  field  because  of  the 
long  thermal-relaxation  time  of  the  cell  in  high  fields.  The  temperature-dependence  was 
fitted  to  the  polynomial  form  P  =  ao  +  a^T  +  02^^  -|-  asT^  in  three  temperature  regions, 
r  <10  mK,  10<  T  <  23  mK,  and  T  >23  mK.  In  this  work  the  region  of  T  <10  mK  was  the 
most  important  because  most  of  the  pressure  change  during  a  sample  precooling  occurred 


^The  temperature-dependence  of  a  filled  cell  could  not  be  measured  because  the  sample  has  its  own 
temperature-dependence. 
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Table  3.1.  Fitting  coefficients  of  the  temperature-dependence  of  the  experimental  cell.  The 
temperature-dependence  has  the  polynomial  form  P  =  ap  +  aiT  +  a2T^  +  asT^-  


ao 

ai 

a-2 

as 

T  <10  mK 

142.689 

-5.15486 

-2.12073 

0.121375 

10  <  T  <  23  mK 

170.696 

-25.6000 

0.971897 

-1.14433x10-'^ 

T  >  23  mK 

-66.4553 

1.01595 

-2.13444x10-'' 

0 

at  these  temperatures.  The  small  temperature-dependence  below  1  mK  was  not  significant 
during  the  demagnetization  of  the  sample. 

3.4    Temperature-Regulation  Method 

The  temperature  of  the  mixing  chamber  was  monitored  and  regulated  using  the  ar- 
rangement shown  in  Fig.  3.9.  The  components  involved  are  a  Linear  Research  110  (LR- 
110)  resistance  bridge,  a  Speer  400  resistance  thermometer  (R4)  located  on  the  mixing 
chamber,  an  Oxford  1  kQ  heater  and  a  LR-130  temperature  controller.  The  value  of  R4 
was  measured  by  the  LR-110  resistance  bridge  whose  DC  error  output  was  used  as  feed 
back  to  the  temperature  controller.  The  LR-130  processed  the  DC  signal  and  delivered  an 
appropriate  amount  of  current  to  the  heater  to  stabilize  the  value  of  R4.  For  this  cryostat 
it  was  found  that  the  best  settings  for  the  LR-110  resistance  bridge  were:  the  excitation 
voltage  Vexc  =  0.3  mVrms]  fine=100%,  and  the  time  constant  r  =  10  seconds.  For  the 
LR-130  temperature  controller,  the  settings  were:  DC  gain  =  30;  variable  integration  time 
constant  Ti  =  30  seconds;  high  frequency  roll-off  time  constant  T2  =  0.03  seconds;  variable 
derivative  time  constant  T3  =  10  seconds;  and  output  resistance  Rout  =  1  kfi  with  the 
monitor  port  shorted  to  bypass  the  1  kfi  current-sensing  resistor. 

3.5    Forming  the  Samples 

All  three  samples,  20.10,  19.65  and  19.61  cm^/mole,  used  the  same  ^He  gas  which  had 
approximately  14  ppm  of  ^He  impurity.  The  impurity  value  was  measured  using  a  Veeco 
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Figure  3.9.  Schematic  drawing  of  the  temperature-control  system. 
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MS-17  leak  detector'''  which  had  been  calibrated  with  a  sample  of  ^He  and  ^He  gas  of  known 
ratio.  The  blocked-capillary  technique  was  used  in  forming  the  samples.  To  grow  a  sample 
at  a  desired  density,  one  must  first  determine  at  what  pressure  the  sample  would  begin 
to  freeze.  For  example,  if  the  desired  molar  volume  was  19.61  cm^/mole  then  the  sample 
would  begin  to  freeze  near  14  MPa  (w  2000  psi)  indicated  by  point  b'  in  Fig.  3.10.  Because 
the  temperature  was  lower  near  the  1  K  pot,  a  solid  block  was  likely  to  form  in  the  fill 
line  before  the  bulk  sample  could  reach  the  melting  curve.  If  this  happened,  one  could  not 
supply  gas  to  the  cell  to  maintain  the  necessary  pressure  and  the  sample  would  reach  the 
melting  curve  at  a  lower  pressure,  thus  yielding  a  lower  density.®  It  was  therefore  necessary 
to  load  the  cell  with  a  pressure  slightly  higher  than  the  targeted  melting  pressure.  The 
loading  pressure  of  our  hep  samples  was  about  14.1  MPa  (2050  psi).  Since  the  solidification 
temperature  of  high-density  solid  ^He  near  the  bcc-hcp  phase  boundary  is  approximately 
3.2  K,  to  form  a  high-quality  sample  requires  good  control  of  the  temperature  in  this  range. 
Below  is  the  description  of  how  high  stability  in  temperature  above  1  K  was  achieved  in 
this  cryostat. 

First,  the  ^He/^He  mixture  was  removed  from  the  dilution  refrigerator.^  Second,  the 
1  K  pot  was  filled  and  pumped  and  used  as  the  cooling  source  in  the  cryostat.  When  cooling 
the  pot  to  1  K,  it  is  pumped  on  by  the  Kinney  pump  through  the  wide  open  #45  valve, 
shown  in  Fig.  2.8.  In  order  to  maintain  higher  temperatures,  the  pumping  rate  must  be 
controlled  using  valve  #47  and  #50,  with  valve  #48  used  as  a  throttle.  For  example,  to 

^Veeco  Instruments,  Inc.,  Terminal  Drive,  Plainview,  NY  11803,  phone:  (516)  349-8300 
*On  this  cryostat,  one  can  melt  the  solid  block  in  the  fill  line  by  applying  60  mW  (10  mA)  to  the  still 
heater. 

^Especially  the  ^He  gas  left  inside  the  refrigerator  would  condense  in  the  heat  exchangers  and  drip  into 
the  mixing  chamber  causing  large  temperature  fluctuations. 
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T  (K) 

\  -  - 

.  ■         I  -  .- 

Figure  3.10.  Molar  volume  and  pressure  versus  temperature  for  high  density  solid  ^He. 
The  line  a-f  represents  a  molar  volume  of  19.61  cm^/mole  with  the  corresponding  pressure 
change  in  the  sample  is  indicated  by  the  path  a'-f . 
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achieve  a  temperature  of  3.2  K  indicated  by  point  b,  b'  in  Fig.  3.10,  the  throttle  valve  was 
opened  wide  enough  so  that  the  average  pressure  of  the  pot  was  about  32  kPa  [44].^° 

Third,  to  cool  other  parts  a  small  amount  of  ^He  gas  was  circulated  in  the  dilution 
refrigerator  using  the  Roots  and  Alcatel  pumps.  The  gas  was  that  remaining  in  the  back 
of  the  Alcatel  pump,  at  a  pressure  of  500  mTorr,  after  all  the  mixture  was  removed  from 
the  refrigerator  and  stored  in  the  storage  tanks. 

When  the  temperature  was  regulated  properly  (see  Section  3.4,  the  refrigerator  was 
left  to  cool  overnight  to  3.4  K,  point  a,  a'  in  Fig.  3.10,  a  temperature  slightly  above  the 
solidification  point  b,  b'  of  the  sample.  For  our  two  hep  samples,  19.61  and  19.65  cm^/mole, 
point  b,  b'  was  about  3.2  K  corresponding  to  a  value  of  1540  Q  of  R4. 

After  the  pressure  in  the  cell  was  set,  the  refrigerator  was  cooled  by  decreasing  the  regu- 
lated temperature  in  0.1  K  (100  Q  on  R4)  steps.  Along  the  melting  curve,  the  temperature 
was  held  fixed  for  5  minutes  at  three  different  values  which  were  about  0.02  K  (AR4  =  20  fi) 
apart,  to  ensure  uniform  sample  density.  Figure  3.11  shows  the  pressure  decrease  during 
the  formation  of  the  19.61  cm^/mole  sample. 

After  the  sample  was  completely  solidified,  it  was  annealed  for  3  days  (just  below  point 
c,  c')  until  the  pressure  decrease  became  negligible.  Figure  3.12  shows  the  pressure  change 
versus  time  during  annealing.  During  the  first  5  hours  the  pressure  increased  by  approx- 
imately 2  kPa  (0.3  psi),  then  began  to  decrease.  For  the  65  hours  of  annealing,  the  total 
pressure  decrease  was  approximately  0.16  MPa  (23  psi).  The  sample  was  then  cooled  slowly 
by  decreasing  the  temperature  in  steps  with  10-minute  waiting  time  per  step.  The  size  of 
the  first  9  steps  was  only  0.1  mK  (1  U),  to  minimize  the  thermal  shock  to  the  sample, 

'"With  these  settings,  the  pressure  of  the  pot  fluctuated  between  17  kPa  (2.8  K)  and  47  kPa  (3.5  K)  with 
a  period  of  approximately  1  hour. 

"The  three  spikes  on  the  melting  curve  are  due  to  a  slight  overshoot  on  the  temperature  controller. 
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Figure  3.11.  Pressure  versus  the  resistance  thermometer  during  the  formation  of  the 
19.61  cm"^/mole  sample.  The  pressure  drop  at  1596  Q  occurs  during  the  annealing  of 
the  sample. 
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Figure  3.12.  Pressure  versus  time  during  annealing  near  the  melting  curve  of  the 
19.61  cm^/mole  sample. 
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following  by  six  0.2  mK  (2  Q)  steps.  Then  the  step-size  was  increased  further  to  0.3,  0.4, 
0.5,  1,  and  2  mK  with  three  steps  each.  The  last  2  mK  step  concluded  the  sample  annealing 
process  for  a  bcc  sample;  however,  a  low-density  hep  solid  requires  another  annealing  which 
occurs  after  the  sample  crosses  the  bcc-hcp  phase  boundary  (see  Fig.  3.10). 

The  bcc-hcp  transition  was  marked  by  a  sharp  drop  in  pressure.  In  our  hep  sample, 
the  drop  was  measured  to  be  approximately  90  kPa  (13  psi),  see  Fig.  3.13.  Additional 
drops  in  pressure  were  achieved  by  changing  the  temperature  of  the  sample  back  and 
forth  below  the  bcc-hcp  transition  which  was  the  procedure  used  for  the  19.65  cm^/mole 
hep  sample.  Alternatively,  for  the  19.61  cm^/mole  sample,  the  same  pressure  drop  was 
achieved  by  lowering  the  temperature  slowly  over  a  period  of  a  week.  The  pressure  drop  at 
a  fixed  temperature  below  the  bcc-hcp  transition  is  shown  in  Fig.  3.14.  The  annealing  data 
indicate  that  the  solid  conversion  from  bcc  to  hep  takes  place  in  a  wide  temperature  range 
below  the  initial  transition.  In  Fig.  3.13,  the  pressure  change  can  be  seen  occurring  with  R4 
as  large  as  1900  Q,  about  150  mK  below  the  bcc-hcp  phase  line.  However,  despite  the  fact 
that  there  was  additional  pressure  drop  during  the  sample  annealing,  the  total  measured 
pressure  decrease  was  only  160  kPa  which  was  200  kPa  smaller  than  what  expected  in 
a  large  bulk  solid  [45].^^  This  phenomena  may  relate  to  the  small  thickness  of  the  bulk 
sample,  which  was  only  0.13  mm.  The  smaller  pressure  decrease  at  the  bcc-hcp  transition 
was  also  reported  by  Okamoto  [27],  whose  sample  thickness  was  0.15  mm.  For  a  sample  of 
molar  volume  of  19.66  cm^/mole,  he  found  a  total  pressure  decrease  of  only  60  kPa  which 
was  250  kPa  smaller  than  the  expected  value  for  a  bulk  sample  [45]. 


The  effect  of  a  non-constant  sample  volume  can  only  accomit  for  approximately  24  kPa  of  the  missing 
pressure  drop. 


Figure  3.13.  Pressure  versus  resistance  thermometer  for  the  19.61  cm^/mole  sample  around 
the  bcc-hcp  phase  line. 
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Figure  3.14.  Pressure  change  at  a  fixed  temperature  during  the  annealing  of  the 
19.61  cm^/mole  sample  below  the  bcc-hcp  phase  line. 


CHAPTER  4 
EXPERIMENTAL  RESULTS 

In  this  chapter,  the  experimental  results  are  presented  with  the  precooling  and  demag- 
netization data  shown  for  a  bcc  sample  of  molar  volume  20.10  cm^/mole  and  two  hep 
samples  with  molar  volumes  of  19.61  and  19.65  cm^/mole.  The  results  are  then  discussed 
in  relation  to  predictions  of  theoretical  models. 

4.1    bcc  Sample 

A  bcc  sample  with  a  molar  volume  of  20.10  cm^/mole  was  used  as  a  test  to  evaluate 
our  experimental  procedure  prior  to  our  measurements  on  hep  samples.  A  bcc  sample  is 
useful  for  this  purpose  because  most  of  its  magnetic  characteristics  have  been  determined 
at  lower  density.  In  addition,  the  pressure  has  been  measured  during  demagnetization  of  a 
bcc  solid  sample  [6]  providing  excellent  data  for  comparison.         •    , ; 

;  V  ■  •    ■    -     '      .  (     -    ■  -  • 

4.1.1    Precooling  Data  '    "  •      -  •       -      •         .  < 

Figure  4.1  shows  the  entire  cooling  process  of  a  ^He  sample.  Our  bcc  sample  was 
precooled  in  a  field  of  2.5  Tesla  over  a  period  of  40  days  to  a  temperature  below  0.55  mK, 
thereby  removing  99%  of  its  entropy,  by  demagnetizing  the  copper  stage.  The  precooling 
process  was  long  due  to  the  large  Kapitza  resistance  [35]  between  the  ^He  solid  and  the  silver 
sinter.  For  magnetic  fields  higher  than  0.2  T,  the  resistance  has  a  l/T^  dependence  which 
gives  rise  to  a  longer  thermal  relaxation  time  of  the  sample  as  the  temperature  decreases. 
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Table  4.1.  Precooling  data  for  v  —  20.10  cm"^/mole  bcc  solid  in  2.5  T  field.  Trundle  is  the 
copper  bundle  temperature  measured  by  the  melting  pressure  thermometer,  is  the  square 
of  the  polarization,  AP  is  the  total  pressure  change  of  the  sample  from  high  temperature, 


p' 

AP  (Pa) 

Waiting  time(hour) 

6.19 

0.091 

-18.6±2.4 

19 

3.2 

0.29 

-43±7.7 

26 

2.55 

0.408 

-131.7±19.1 

115 

1.75 

0.642 

-129.6±7.2 

44 

0.88 

0.952 

-317.1±8.9 

121 

0.55 

0.990 

-284±1.1 

235 

However,  this  large  resistance  is  beneficial  in  preventing  a  large  amount  of  heat  flow  from 
the  cell  into  the  sample  during  the  demagnetization  of  the  sample. 

The  precooling  of  the  sample  was  carried  out  in  steps.  First,  the  sample  was  precooled 
for  107  hours^  along  with  the  copper  bundle  to  a  temperature  near  11  mK.  Second,  the 
copper  bundle  was  demagnetized  to  6  mK  where  the  temperature  was  held  constant  for 
19  hours  allowing  the  sample  to  cool  to  the  new  temperature.  Then  the  temperature 
of  the  copper  bundle  was  reduced  further,  with  this  process  continuing  until  the  sample 
reached  0.55  mK  as  shown  in  Table  4.1.  At  each  temperature,  the  pressure  change  AP  was 
determined  by  fitting  to  the  formula: 


AP=  Po  +  Piexp(-Vr)  (4.1) 

where  Pq  +  Pi  is  the  pressure  at  t  =  0  when  the  bundle  reaches  the  new  temperature,  and 
r  is  the  thermal  time  constant. 

The  pressure  change  in  the  sample  versus  the  square  of  the  polarization  which  is  given 
by  p  =  ta.nh{yhBo/2kBT)  and  the  temperature  are  plotted  in  Figs.  4.2  and  4.3,  respectively. 

'A  typical  precooling  takes  about  72  hours;  however,  in  this  run  our  dilution  refrigerator  was  partially 
blocked  and  the  circulation  rate  was  reduced  to  maintain  an  operational  return-line  pressure. 
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Figure  4.2.  Pressure  change  versus  the  square  of  the  polarization  during  the  bcc  sample 
precooling. 

The  solid  lines  in  both  graphs  were  obtained  from  Okamoto's  work  on  a  bcc  sample  with 
a  molar  volume  of  19.99  cm^/mole  [6]  scaled  to  the  molar  volume  of  our  sample.  While 
our  objective  was  not  to  determine  the  precooling  with  high  accuracy,  our  precooling  data 
agree  well  with  those  of  the  Tokyo  group.  Using  the  data  in  Fig.  4.2  together  with  Eq.  1.35, 
one  can  determine  the  volume  derivatives,  dByj/dv  and  dKp/dv  for  a  bcc  solid.  For  their 
bcc  solid,  Okamoto  et  al.  [6]  found  d9^/dv  =  -40.7  ±  2.2  /iK  mole/cm^  and  dK'^/dv  = 
-11.0  ±  0.7 //K  mole/cm^. 
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Figure  4.3.  Pressure  change  versus  temperature  during  the  bcc  sample  precooling. 
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4.1.2  Demagnetization  Data 

After  precooling  in  2.5  T  to  0.55  mK,  the  sample  was  demagnetized  with  a  ramp-rate 
of -73  mT/hour.  This  rate  was  kept  constant  to  zero  field  in  order  to  avoid  the  spikes  that 
occur  with  changes  in  ramp-rate  (see  Fig.  3.6).  Figure  4.4  shows  the  pressure  curve  of  the 
bcc  sample  during  its  demagnetization,  with  the  field-dependence  of  the  cell  subtracted.  A 
similar  pressure  behavior  was  observed  by  Okamoto  et  al  [6].  The  pressure  remained  flat 
until  0.5  T  and  then  began  to  rise,  indicating  that  the  sample  had  entered  the  high-field 
phase  from  the  paramagnetic  phase.  The  pressure  increased  rapidly  between  the  fields 
of  0.5  T  and  0.02  T  with  a  total  change  of  340  Pa.  A  second  transition,  marked  by  a 
pressure  jump  of  20  Pa,  occurred  at  0.02  T  when  the  sample  entered  the  LFP  from  the 
HFP.  The  observation  of  this  transition  indicates  that  our  bcc  sample  was  cooled  to  a 
temperature  below  30  fiK.  The  pressure  of  our  sample  remained  flat  in  the  LFP  to  zero 
field,  unlike  the  Tokyo  data  [6]  which  exhibited  a  pressure  decrease  of  35  Pa  in  the  LFP. 
The  absence  of  a  pressure  drop  at  low  fields  indicates  that  the  use  of  a  Stycast  seal  instead 
of  an  indium  0-ring  in  our  sample  cell  significantly  reduced  heating  near  the  end  of  the 
sample  demagnetization.  This  knowledge  is  important  in  explaining  the  results  of  our  hep 
samples. 

4.1.3  Discussion  of  bcc  Results 

To  understand  better  the  bcc  pressure  data,  it  is  useful  to  examine  a  few  adiabatic 
demagnetization  paths  on  the  B  -  T  phase  diagram  for  different  entropies.  For  each  of 
these  paths,  the  value  of  the  entropy  was  determined  from  the  initial  temperature  and  field 
at  the  beginning  of  the  demagnetization.  The  B  -  T  phase  diagram  was  constructed  for 
the  bcc  sample  using  the  results  from  two  previous  experiments.  The  LFP-PP,  LFP-HFP 


92 


200 


100 


-I  I  I  I  I  I  I  I  r 


120 


—  low-field 
*d*  phase 


100 


<! 


-100 


LFP      !  HFP 


80 
0.00 


0.05 


B  (T) 


-200 


-300 


-400 


high-field  phase 


-1  I  I  L 


paramagnetic  phase 


J  I  L 


0.0 


0.2 


0.4  0.6 

B  (T) 


0.8 


1.0 


Figure  4.4.  Pressure  change  during  the  bcc  sample  demagnetization  with  nearly  100% 
polarization. 
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and  low-field  PP-HFP  boundaries  were  scaled  from  the  work  of  Xia  et  al.  [13]  at  melting 
density  using  scaling  factors  of  v^^  for  temperature  and  t;^^-^  for  field  [12].  The  high-field 
PP-HFP  boundary  was  scaled  from  work  of  Fukuyama  et  al.  [46]  in  which  the  magnetic 
phase  diagram  of  the  bcc  solid  at  a  molar  volume  of  22.69  cm^/mole  was  determined  by 
pressure  measurements.  The  scaling  factor  used  for  both  the  temperature  and  the  fields 
was  u^®.  (,  J 

Adiabatic  demagnetization  curves  were  generated  and  scaled  to  our  sample  using  the 
constant-field  entropy  data  at  melting  density  published  by  Ni  et  al.  [29].  The  scaling  factors 
were  v^^  and  u^'*'^  for  the  temperatures  and  the  fields,  respectively.  Figure  4.5  shows  three 
predicted  paths  on  a.  B  —  T  phase  diagram  for  adiabatic  sample  demagnetizations  of  the  bcc 
solid  with  S/R\n2  =  0.1,  0.3,  and  0.7.  The  entropy  values  can  be  determined  by  substituting 
in  the  initial  values  of  the  temperature  and  applied  field  into  the  equation: 

5  =  i2(ln(2coshx)  -  a;tanhx),  (4.2) 

where  R  is  the  gas  constant;  a;  is  a  dimensionless  parameter  given  by  fhBo/2kBT  with  Bq 
the  applied  field. ^ 

For  S  =  0.1ir:in2  the  sample  enters  the  HFP  at  about  0.3  T,  with  the  cooling  rate 
decreasing  with  decreasing  field.  The  HFP-LFP  boundary  is  reached  near  0.03  T  and 
21  fiK.  Then,  because  of  the  lower  entropy  in  the  LFP,  the  sample  warms  slightly  along 
this  boundary  to  a  temperature  of  26  /xK  where  it  enters  the  LFP  and  remains  to  zero  field. 
The  demagnetization  path  of  our  bcc  sample  was  similar  to  that  for  S  =  0.1Mn2,  except 
that  our  sample  entered  the  HFP  at  slightly  higher  field  because  its  remaining  entropy 

^The  above  equation  is  only  valid  in  high  magnetic  fields  where  the  Zeeman  energy  is  dominant. 
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was  almost  0.  There  is  also  a  slight  disagreement  in  the  value  of  Bd  (HFP-LFP).  In  our 
measurement,  we  found  Bd  to  be  about  0.02  T  yielding  a  scaling  factor  of  16.7  which  is 
higher  than  the  14.5±0.3  obtained  in  Fukuyama's  work  [12].  For  their  bcc  sample  with 
a  molar  volume  of  19.99  cm^/mole,  Okamoto  et  al.  [6]  reported  Bd  to  be  0.027  T  which 
yielded  a  scaling  factor  of  about  14.7.  However,  because  the  applied  magnetic  field  was 
trapped  when  the  indium  0-ring  became  superconducting  {Be  =  0.029  T),  Okamoto  and 
co-workers  were  forced  to  measure  the  trapped  field  in  a  dummy  cell  and  then  assumed 
that  the  same  field  was  trapped  in  the  real  cell  in  order  to  deduce  the  value  of  the  field 
at  the  sample.  In  our  experiment,  the  trapped-field  problem  was  totally  avoided  because 
our  cell  has  no  indium  0-ring  and  the  value  of  the  field  was  measured  directly.  Therefore, 
the  measured  value  for  Bd  in  our  sample  is  more  reliable  than  that  of  Okamoto  et  al. 
who  might  have  over-estimated  the  value  of  the  trapped  field  in  their  cell.  The  question  of 
the  scaling  factor  for  Bd{v)  bears  further  investigation  because  of  its  implication  for  the 
volume-dependence  of  various  exchange  parameters. 

The  second  demagnetization  curve  is  for  S  =  0.3/21n2.  In  the  HFP,  the  behavior  is 
almost  like  that  for  S  =  0.1i?ln2;  however,  the  sample  never  enters  the  LFP  fully  but  follows 
the  HFP-LFP  and  then  the  PP-LFP  phase  lines  to  zero  field.  In  other  words,  the  sample  is 
only  partially  ordered  at  low  fields.  The  third  curve  represents  the  demagnetization  curve 
of  a  sample  with  S  =  0.7i21n2.  The  sample  neither  enters  the  HFP  nor  the  LFP,  remains 
in  the  PP  even  at  zero  field. 

The  predicted  path  for  S  =  0.1i?ln2  is  qualitatively  in  agreement  with  our  bcc  exper- 
imental data.  Both  results  clearly  show  two  magnetic  ordering  transitions  (PP-HFP  and 
HFP-LFP)  and  the  observation  of  the  HFP-LFP  transition  indicates  that  the  sample  was 
cooled  below  30  /xK.  The  predicted  paths  also  show  that  magnetic  ordering  can  be  achieved 
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only  at  a  certain  value  of  entropy,  and  different  initial  conditions  will  yield  different  results. 
This  knowledge  is  extremely  useful  in  examining  the  pressure  data  of  hep  samples. 

4.2    hep  Sample 

In  two  separate  coolings  from  room  temperature  in  March  1995  and  October  1995,  pres- 
sure measurements  were  made  on  two  hep  samples  with  molar  volumes  of  19.65  cm"^/mole 
and  19.61  em^/mole,  respectively.  This  section  is  divided  into  two  parts,  one  for  the  pre- 
cooling  data  and  the  other  for  the  demagnetization  data.  There  were  a  total  of  seven 
precoolings  and  demagnetizations,  three  for  the  first  sample  and  four  for  the  second,  with 
conditions  as  described  in  the  following  text. 

4.2.1    Preeooling  Data 

In  the  first  preeooling,  the  hep  sample  (i;  19.65  em^/mole)  was  cooled  to  a  temperature 
of  0.61  mK  while  in  a  field  of  2.5  T.  The  duration  of  the  entire  process  was  38  days  and  at 
the  end  the  spins  were  almost  fully  polarized  with  about  98%  of  the  spin  entropy  removed 
from  the  sample. 

As  in  the  preeooling  of  the  bee  sample,  the  hep  preeooling  was  carried  out  by  slowly 
stepping  down  the  temperature  of  the  copper  bundle.  In  the  first  step,  the  sample  was 
precooled  by  the  dilution  refrigerator,  along  with  the  copper  bundle,  to  11.2  mK  in  84 
hours.  The  temperature  of  the  bundle  was  then  lowered  further  by  demagnetization  as 
shown  in  Table  4.2.  As  discussed  in  the  previous  section  for  the  bee  sample,  the  pressure 
change  of  the  hep  sample  for  each  preeooling  steps  was  fitted  to  an  exponential  form. 
Figure  4.6  shows  the  plot  of  AP  versus  of  the  sample  during  three  different  precoolings. 
Here,      was  determined  by  the  applied  field  Bq  and  the  preeooling  temperature  T  using 
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Table  4.2.  First  precooling  data  for  v  =  19.65  cm^/mole  hep  solid  in  2.5  T  field.  Trundle 
is  the  temperature  of  copper  bundle  measured  by  the  melting  pressure  thermometer,  is 
the  square  of  the  polarization  of  the  spins,  AP  is  the  total  pressure  change  of  the  sample 
from  high  temperature,  and  the  waiting  time  is  the  amount  of  time  during  which  Tbundle  is 


Tbundle{mK) 

p' 

AP  (Pa) 

Waiting  time(hour) 

6.06 

0.095 

13.1±1.6 

24 

3.38 

0.266 

24.6±4.6 

85 

2.45 

0.431 

55.5±13.5 

71 

1.78 

0.631 

36.1±19.3 

22 

0.954 

0.932 

147.5±10.3 

142 

0.614 

0.993 

119.6±0.8 

92 

Table  4.3.  Second  precooling  data  for  v  —  19.65  cm"^/mole  hep  solid  in  2.5  T  field. 


Tbundle  i^^) 

AP  (Pa) 

Waiting  time(hour) 

6.73 

0.078 

10.8±9.5 

19 

3.61 

0.239 

18.3±4.6 

64 

2.79 

0.358 

22.1±21.3 

21 

the  equation  p  —  ta,nh{jhBo/2kBT).  The  circles  represent  the  first  precooling  data  and  the 
solid  line  is  the  fitted  curve  extrapolated  from  the  Tokyo  data  [6]  for  this  sample. 

In  the  second  and  third  precoolings,  the  sample  was  cooled  to  a  temperature  near 
3  mK  in  a  2.5-T  field.  The  amount  of  entropy  remaining  in  the  system  was  approximately 
0.75i?ln2  with  a  polarization  of  50%.  The  reason  for  leaving  the  higher  entropy  in  the 
sample  is  discussed  in  Section  4.2.2.  The  two  precooling  sequences  are  given  in  Tables  4.3 
and  4.4.  '  ' '    "  .  '  i  '  ! 

For  the  second  hep  sample  with  t;  =  19.61  cm^/mole,  there  were  four  precoolings  (the 
fourth,  fifth,  sixth,  and  seventh  of  the  precoolings  of  hep  samples).  In  the  fourth  precooling. 


Table  4.4.  Third  precooling  data  for  v  =  19.65  em^/mole  hep  solid  in  2.5  T  field. 


Tbundle  (mK) 

P' 

AP  (Pa) 

Waiting  time(hour) 

7.02 

0.072 

10.0±9.6 

19 

3.90 

0.209 

45.6ib4.7 

96 

3.00 

0.321 

35.1±25.1 

19 

160 


Figure  4.6.  Pressure  change  versus  the  square  of  the  polarization  during  the  sample 
cooling  of  the  first  hep  sample. 
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Table  4.5.  Fourth  precooling  data  for  v  =  19.61cm'^/mole  hep  solid  in  2.5  T  field. 


T6u„d;e(niK) 

AF  (Pa) 

Waiting  time{hour) 

6.13 

0.093 

12.9±1.1 

46 

3.99 

0.201 

39.1±4.0 

90 

Table  4.6.  Fifth  precooling  data  for  v  =  19.61cm'^/mole  hep  solid  in  2.5  T  field. 


Tbundle{mK) 

p' 

AF  (Pa) 

Waiting  time(hour) 

6.05 

0.095 

13.1±1.2 

65 

3.48 

0.254 

23.9±1.6 

116 

1.77 

0.635 

58.1±0.8 

325 

the  sample  was  cooled  to  4  mK  with  a  remaining  entropy  of  0.85i21n2.  In  the  fifth  precooling, 
the  sample  had  a  final  temperature  of  1.77  mK  corresponding  to  a  remaining  entropy  of 
0.46i?ln2.  The  sixth  and  seventh  precoolings  had  the  same  final  temperature  of  1.20  mK 
corresponding  to  0.23i?ln2  entropy.  The  seventh  precooling  was  a  repetition  of  the  sixth 
because  of  a  large  heat  leak  encountered  during  the  sixth  demagnetization  of  the  hep  sample. 
Tables  4.5,  4.6,  and  4.7  list  the  sequences  for  the  fourth,  fifth,  and  sixth  precoolings  of  the 
hep  sample,  respectively.  Figure  4.7  shows  a  plot  of  AF  versus  for  the  second  sample 
with  three  different  symbols  representing  the  three  different  precoolings. 

4.2.2    Demagnetization  Data 

In  this  section,  the  first  observation  of  magnetic  ordering  of  a  hep  sample  is  presented. 
Our  results  indicate  that  the  pressure  of  the  sample  behaves  differently  depending  on  the 
amount  of  entropy  removed  during  the  precooling.  By  comparing  our  data  for  different 
entropies  to  calculations  for  the  high-temperature  expansion  of  the  multiple-exchange  and 

Table  4.7.  Sixth  precooling  data  for  v  =  19.61em^/mole  hep  solid  in  2.5  T  field. 


76und;e(mK) 

P' 

AF  (Pa) 

Waiting  time(hour) 

6.36 

0.087 

ll.lil.O 

41 

1.74 

0.645 

91.1±23.5 

114 

1.20 

0.851 

83.1±0.7 

472 

Figure  4.7.  Pressure  change  versus  the  square  of  the  polarization  during  the  sample  pre- 
cooling  of  the  second  hep  sample. 
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Figure  4.8.  Demagnetization  data  for  the  v  =  19.65  cm^/mole  sample  with  entropy  of 
0.02i?ln2. 

Heisenberg  models  discussed  in  Section  4.2.3,  we  have  strong  evidence  that  our  low-entropy 
hep  sample  was  magnetically  ordered  at  low  temperature.  The  demagnetization  data  for 
the  first  sample,  v  =  19.65  cm-^/mole,  is  discussed  first  followed  by  a  discussion  of  the 
results  from  the  second  sample,  v  =  19.61  cm^/mole. 

The  first  demagnetization  of  the  hep  samples  (v  =  19.65  cm^/mole)  started  with  98% 
entropy  removal.  The  change  in  pressure  versus  field  during  the  demagnetization  is  shown 
in  Fig.  4.8.  The  pressure  increased  by  approximately  6  Pa  over  the  field  range  from  1.8  T 
to  0.2  T  and  then  began  to  drop,  with  a  total  amount  of  10  Pa  between  0.2  T  and  0  T.  A 
similar  pressure  decrease  in  low  fields  was  also  observed  by  the  Tokyo  group  [6];  however, 
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they  attributed  this  pressure  drop  to  heating  of  their  hep  samples.  In  our  case  no  heating 
occurs  at  low  fields  as  demonstrated  for  the  bcc  solid. 

In  a  previous  publication  [47],  we  raised  the  possibility  that  the  decrease  in  pressure  was 
the  signature  of  ordering  of  the  hep  solid.  To  demonstrate  that  this  is  the  case,  we  carried 
out  two  tasks.  First,  we  made  pressure  measurements  on  samples  with  high  entropies  so  that 
the  samples  could  not  magnetically  order  near  the  end  of  a  demagnetization.  These  results 
showed  the  behavior  of  unordered  samples  during  a  sample  demagnetization  and  allowed  us 
to  make  a  comparison  with  the  results  for  a  low-entropy  sample  in  which  magnetic  ordering 
occurred.  Second,  we  used  theoretical  models  for  comparison  with  our  pressure  measure- 
ments on  high-entropy  samples.  For  the  second  task,  we  employed  the  high-temperature 
expansion  of  the  multiple-exchange  and  Heisenberg  models,  described  in  Section  4.2.3.  The 
pressure  measurements  of  high-entropy  samples  are  discussed  below. 

The  second  precooling  ended  when  the  sample  (u  =  19.65  cm^/mole)  reached  a  temper- 
ature of  approximately  2.8  mK,  corresponding  to  a  remaining  entropy  of  0.72i21n2.  During 
the  sample  demagnetization,  shown  in  Fig.  4.9,  the  pressure  increased  by  approximately 
11  Pa  between  1.8  and  0.1  T,  a  behavior  quite  similar  to  that  of  the  first  precooling  in 
which  5  «  0.  However,  below  0.1  T,  the  pressure  increased  rapidly  with  a  total  change  of 
26  Pa  in  going  to  zero  field  (from  ~  34  Pa  at  0.1  T  to  ~  60  Pa  at  0  T).  The  increase  in 
pressure  was  quite  a  contrast  to  the  behavior  observed  in  the  first  sample  demagnetization. 
The  implication  of  this  different  behavior  will  be  discussed  in  Section  4.2.3. 

Our  third  hcp-sample  demagnetization  was  with  0.75/?ln2  remaining  entropy  and  we 
observed  a  similar  pressure  increase.  The  only  major  diflFerence  between  the  second  and  the 
third  demagnetization  was  the  temperature  of  the  copper  bundle,  which  was  kept  at  0.5  mK 
near  the  end  of  the  third  demagnetization  instead  of  0.1  mK  for  the  second  demagnetization. 
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Figure  4.9.  Demagnetization  data  for  the  v  =  19.65  cm^/mole  sample  with  entropy  of 
0.72i?ln2. 
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Figure  4.10.  Demagnetization  data  for  the  v  =  19.61  cm^/mole  sample  with  entropy  of 
0.85flln2. 

The  final  bundle  temperature  in  the  first  demagnetization  was  also  0.5  mK  and  we  wanted  to 
be  sure  that  the  pressure  increase  that  we  first  encountered  in  the  second  demagnetization 
was  not  related  to  the  temperature  of  the  copper  bundle. 

The  fourth  sample  demagnetization  (u  =  19.61  cm^/mole)  was  with  5  =  0.85i?ln2.  As 
shown  in  Fig.  4.10,  the  pressure  of  the  sample  rose  approximately  20  Pa  towards  the  end 
of  the  demagnetization.  As  will  be  discussed  in  Section  4.2.3,  this  experimental  result  also 
agreed  with  the  S  =  0.85i21n2  curve  in  the  exchange-model  calculation. 

The  fifth  sample  demagnetization  was  carried  out  with  5  =  0.46i?ln2  and  yielded  no 
significant  pressure  change  down  to  zero  field  as  shown  in  Fig.  4.11.   This  is  a  strong 
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Figure  4.11.  Demagnetization  data  of  the  v  —  19.61  cm^/mole  sample  with  entropy  of 
0.46i?ln2. 

indication  that  the  sample  was  partially  ordered  since  the  pressure  change  was  between 
those  of  the  0  and  0.75i21n2  results. 

The  sixth  sample  demagnetization  had  an  initial  entropy  of  0.23i21n2  but  was  ruined  by 
a  large  heat  leak  of  approximately  100  nW  to  the  copper  bundle.  Despite  an  attempt  to  cool 
the  bundle  by  demagnetization,  the  bundle  warmed  up  from  1  mK  to  about  5  mK  during 
the  34  hour  sample  demagnetization.  Two  probable  causes  for  the  large  heat  leak  were 
the  power-supply- ramp  unit  used  in  the  demagnetization  of  the  sample  or  the  defective 
lock-in  amplifier  of  the  melting  pressure  thermometer.  The  thermometer  was  unusually 
noisy  during  the  precooling  of  the  sample.  Figure  4.12  shows  the  sixth  demagnetization 
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Figure  4.12.  Demagnetization  data  for  the  t;  =  19.61  cm^/mole  sample  with  large  heat 
leak. 

curve  together  with  the  fourth  and  fifth  demagnetization  data.  Because  of  its  initial  lower 

;     >  •     >  «         .     t  .■     ,  J 

«  '■.•■>•"'.'-     ^  . 

entropy,  the  sample  pressure  in  the  sixth  demagnetization  was  higher  than  those  of  the 
fourth  and  fifth.  However,  as  more  heat  flowed  into  the  sample,  the  entropy  increased 
and  the  pressure  dropped  with  decreasing  field.  At  B  =  0.45  T,  the  sample  entropy  had 
increased  to  about  0.46i?ln2  and  to  0.85i?ln2  at  0.15  T,  becoming  almost  i?ln2  at  zero  field. 

The  seventh  sample  demagnetization,  as  shown  in  Fig.  4.13,  was  a  repetition  of  the 
sixth  with  S  =  0.23i?ln2.  The  overall  pressure  change  during  the  demagnetization  is  sim- 
ilar to  that  of  the  first  demagnetization  with  S  =  0.02i?ln2.  In  both  cases,  the  pressure 
measurements  show  a  slight  increase  from  high  fields  down  to  approximately  0.2  T,  then  a 
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Figure  4.13.  Demagnetization  data  of  the  v  =  19.61  cm^/mole  sample  with  entropy  of 
0.23i?ln2. 

decrease  down  to  zero  applied  field.  The  only  difference  between  the  two  cases  was  the  size 
of  the  pressure  decrease  in  low  fields.  For  the  seventh  sample  demagnetization,  the  total 
pressure  drop  was  approximately  3  Pa  between  0.2  T  and  0  T  compared  to  10  Pa  observed 
in  the  first  demagnetization.  As  will  be  presented  in  Section  4.2.3,  this  experimental  result 
disagrees  with  the  high- temperature  calculations  for  the  multiple-exchange  and  Heisenberg 
models. 

The  demagnetization  data  clearly  show  that  the  pressure  of  the  hep  samples  behaves 
differently  depending  on  the  amount  of  entropy  remaining  in  the  sample.  High-temperature 
expansions  of  the  multiple-exchange  and  Heisenberg  models,  which  are  discussed  next,  give 


108 


strong  support  that  pressure  increase  at  low  fields  of  high-entropy  samples  is  the  result 
of  adiabatic  demagnetization  of  a  disordered  ferromagnet.  For  the  low-entropy  sample, 
the  best  explanation  for  the  pressure  decrease  at  low  field  is  the  ferromagnetic  ordering 
transition  of  the  hep  solid. 

4.2.3    Comparison  of  Results  with  Model  Calculations 

Two  theoretical  models,  multiple-exchange  and  Heisenberg,  were  used  in  simulation 
of  the  sample  demagnetizations.  The  main  purpose  of  the  simulation  was  to  predict  the 
pressure  change  in  an  hep  sample  during  a  direct  demagnetization  for  comparison  with  the 
experimental  results.  These  are  high-temperature  expansions  valid  only  for  T  >  Tc,  the 
ordering  temperature. 

Exchange  Model 

The  calculation  is  based  on  the  expression  for  the  free  energy  provided  by  Roger  et  al. 
[23]  of  their  high-temperature  expansion  of  the  two-,  three-,  and  four-particle  exchange 
model  (for  more  detail  see  Section  1.2.2).  The  expansion,  for  a  sample  located  in  an 
arbitrary  field,  was  carried  to  the  fourth  order  of  l/ksT.  The  expression  for  the  free  energy 
for  an  hep  solid  is  given  by 

(4.3) 

where  R  is  the  gas  constant;  a;  is  a  dimensionless  parameter  given  by  jhBQ/2kBT,  with  Bq 
the  applied  field;  p  (=  tanh  x)  is  the  polarization;  A;b^^  =  3Ji  +  SJ2  +  9K'^^  with  Ji  and  J2 
the  in-plane  nearest-neighbor  and  out-of-plane  nearest-neighbor  contributions  of  the  two- 
particle  and  three-particle  exchanges,  respectively,  and  K'^^  the  four-particle  out-of-basal 
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plane  exchange;  and  the  an(p^)  are  polynomials  of  order  2n  -  1  in  whose  coefficients  are 
shown  in  Table  4.8.  By  taking  the  derivative  with  respect  to  the  temperature  at  constant 
field  and  molar  volume,  —dF/dT\Bo,v,  one  obtains  an  expression  for  the  entropy 

S   -    R  {\n{2  cosh x)  -  xp  - -—^^—^(kBe^p  +  SK'  p^)  (4.4) 
L  Kfijcosh  X  ^  ' 

2xp    J^anjp^)  (   1  V 

^,^cK2  ^  on 


cosh^x        2"n!  \kBT 

n=z 

-  (i-/)E 


n=2 


X        dan{p^)       ,  .      ,  2s 

+  (n  -  l)a„(p  ) 


cosh^(x)  dp 


2"n!(A;Br)"  j  ' 


where  dan{p^)/dx  is  the  partial  derivative  of  the  polynomials  a„(p^)  with  respect  to  x. 

The  expression  for  the  pressure  can  be  obtained  from  the  free  energy  by  taking  its 
derivative  with  respect  to  volume, 


P  =  R 


dv  2  ^      dv       4^    ^      ^'^^      2"n!  UbT/ 


(4.5) 


where  the  dan{p^)/dv  is  the  partial  derivative  of  the  polynomials  an(p^)  with  respect  to 
the  molar  volume  v. 

A  computer  program  was  written  to  calculate  the  pressure  change  during  an  adiabatic 
demagnetization  of  a  hep  sample.  The  program  first  calculates  the  entropy  for  given  initial 
temperature  and  field.  The  field  was  then  reduced  in  steps  to  simulate  the  demagnetization. 
At  each  new  value  of  the  field  a  root-finding  routine,  "Newton-Raphson  and  Bisection"  [48], 
was  applied  to  Eq.  4.4  with  the  entropy  held  constant,  in  order  to  determine  the  new  tem- 
perature T.  This  new  value  was  substituted  into  Eq.  4.5  to  calculate  the  pressure  change. 
Figure  4.14  shows  the  results  of  the  calculation  at  different  entropies  for  an  hep  sample  of 
molar  volume  of  19.65  cm^/mole.  Since  the  model  breaks  down  at  low  temperatures  (fields). 
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Table  4.8.  Polynomial  coefficients,  obtained  from  Ref.  [23],  for  Aa^-y(p^)  of  the  exchange 
model  forn  =  2  to  4.  Aa/3^(p^)  are  related  to  the  polynomial  a„{p^)  of  the  high-temperature 
expansion  by  a„(p^)  =  ^  JfJf ''^'^^a07(p^)-  For  a  given  n  and  each  combination  of 
Jf  Jf/f"^,  the  coefficients  are  listed  in  columns  with  increasing  power  of  p^  up  to  (p^)^"~^ 
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Figure  4.14.  Calculated  pressure  change  relative  to  the  unpolarized  state  using  the  high- 
temperature  expansion  of  the  multiple-exchange  model  during  a  sample  demagnetization 
at  different  initial  entropies. 
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the  calculation  in  this  range  is  excluded  from  the  graph.  Except  for  only  a  small  pressure 
increase  for  the  lowest  entropy,  all  curves  show  large  upward  turns  of  the  pressure  at  low 
fields.  Our  experimental  data  with  S  =  0.72i?ln2  and  0.85i21n2  agreed  qualitatively  with 
the  calculation;  however,  the  pressure  of  a  lower-entropy  samples  did  not  behave  according 
to  the  theoretical  prediction.  We  interpret  this  to  mean  that  the  condition  T  >  Tc  required 
for  validity  of  the  expansion  is  not  met  because  Tc  has  been  reached. 

In  the  simulation  we  chose  Ji  =  J2  as  assumed  by  Roger  et  al.  [23]  and  then  used 
the  precooling  pressure  data  of  the  Tokyo  group  [6]  to  determine  the  value  of  K'^^.  For 
their  hep  sample  of  u  =  19.66  cm"^/mole,  they  found  dOyj/dv  =  33.6  mole/cm^  and 
\/kB{dK[Jdv)  =  -3.8  /iK  mole/cm^.  Since  ksOw  =  3Ji  +  3J2  +  ^K'    we  have 


The  above  result  indicates  that  the  value  for  K'^^  from  the  pressure  measurement  is  ap- 
proximately 5/3  times  larger  than  the  estimation  for  K'^^  by  Roger  et  al.  [23]  who  used  the 
susceptibility  data  at  v  =  19.4cm3/mole  [2].  Other  susceptibility  data  [49]  and  calculations^ 
yields  Ji  =  -3.91/^^,  and  -4.82A'^,  respectively  for  v  =  19.4  cm^/mole.  We  believe  that 
the  pressure  data  yields  a  better  estimation  since  they  were  obtained  at  high  temperature 
(during  the  precooling  of  the  sample)  where  the  expansion  is  valid.  The  use  of  the  sus- 
ceptibility data  for  estimation  of  the  exchange  parameters  is  less  preferable  because  one  is 
forced  to  make  the  fitting  at  a  very  low  temperature  where  the  high-temperature  expansion 
may  not  be  valid.  ,  v     .  '1: 

'D.  Ceperley,  private  communication. 
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Figure  4.15.  Pressure  change  (solid  dots)  relative  to  the  unpolarized  state  during  a  sample 
demagnetization  at  different  initial  entropies  and  calculation  from  the  high-temperature 
expansion  of  the  multiple-exchange  model  (lines)  for  v  =  19.65  cm^/mole. 

Comparing  this  experimental  result  with  that  for  the  S  =  0.72Mn2  high-temperature 
expansion  in  Fig.  4.15,  we  find  that  both  of  them  have  increasing  pressures  for  decreasing 
field  at  low  fields.  The  magnitude  of  the  pressure  rise  and  the  field  value  at  which  this  rise 
occurs  are  also  similar.  These  calculated  results  provide  strong  evidence  that  the  observed 
pressure  increase  is  the  result  of  adiabatic  demagnetization  of  a  disordered  ferromagnet. 
For  the  low-entropy  sample  with  S  =  0.02i21n2,  the  multiple-exchange  model  also  predicts 
a  small  pressure  increase;  however,  that  was  not  observed  for  that  case.  The  propable  ex- 
planation is  that  the  decrease  in  pressure  is  because  of  the  ferromagnetic  ordering  signature 
of  the  hep  solid  ^He.  The  high-temperature  series  expansion  would  not  be  valid  there. 
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Figure  4.16.  Pressure  change  (solid  dots)  during  a  sample  demagnetization  at  different 
initial  entropies  and  calculations  from  the  high-temperature  expansion  of  the  exchange 
model  (lines)  for  t;  =  19.61  cm^/mole. 

Figure  4.16  shows  the  plot  of  the  calculated  pressure  together  with  experimental  data 
for  the  second  hep  sample  v  =  19.61  cm^/mole.  The  data  for  S  =  0.85Mn2  exhibited 
an  increase  in  pressure  at  low  fields,  qualitatively  in  agreement  with  the  calculated  result 
although  not  quite  as  large.  The  disagreement  in  the  pressure  magnitude  may  be  due 
to  the  fact  that  the  high-temperature  expansion  of  the  model  has  been  carried  only  to  the 
fourth  order  of  l/ksT  or  because  of  incorrect  estimation  in  the  values  of  the  three  exchange 
parameters  Ji,  J2,  and  K'^^.  Our  simulation  using  the  high-temperature  expansion  up  to 
10th  order  of  l/ksT  of  the  Heisenberg  model  (discussed  in  the  next  section)  yielded  a 
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smaller  change  in  pressure  at  low  fields  than  that  obtained  from  the  multiple-exchange 
model  for  the  same  entropy  value.  ^  ,  -  "  , 

For  the  lower-entropy  data  of  5  =  0.46  and  0.23i21n2,  there  was  large  disagreement  with 
the  calculated  results.  The  S  =  0.46i?ln2  data  had  only  a  small  pressure  increase  of  approx- 
imately 5  Pa  at  low  fields;  however,  the  calculated  pressure  change  for  the  same  entropy 
yielded  a  positive  change  in  AP  of  more  than  40  Pa.  For  S  =  0.23/21n2,  the  experimental 
data  had  a  3  Pa  pressure  decrease  near  zero  field  while  the  exchange  model  predicted  a 
positive  pressure  change  of  more  than  40  Pa.  These  large  disagreements  indicated  that 
there  was  magnetic  ordering  of  the  nuclear  spins. 

Because  of  large  uncertainty  in  the  pressure  during  the  sample  precoolings  (Figs.  4.6 
and  4.7),  the  experimental  data  have  been  shifted  so  that  they  agree  with  the  calculation 
at  high  fields  where  the  calculated  pressure  change  is  insignificant.  In  both  figures,  the 
measurements  of  high-entropy  samples  qualitatively  agree  with  the  calculated  results.  For 
lower-entropy  samples,  the  calculation  fails  to  predict  the  pressure  change  at  low  fields.  The 
discrepancy  becomes  larger  with  decreasing  entropy,  giving  strong  evidence  that  nuclear 
spins  in  the  sample  are  more  strongly  ordered  at  low  entropy. 

Heisenberg  One-Exchange  Model 

To  check  for  the  validity  of  our  multiple-exchange  model  simulation  and  to  extend  our 
calculation  to  lower  temperature  (field),  we  employed  the  high-temperature  expansion  of 
the  Heisenberg  model  up  to  the  tenth  order  in  J/ksT  where  J  is  the  interaction  of  two 
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nearest  neighbors  [50].  The  free  energy  for  this  model  is  given  by"* 


F^-RT 


ln(2coshx)  +  5:  Yl 


n>l 0<j<n 


(4.7) 


where  x  =  hjBo/2kBT  and  p  =  tanha;  are  the  same  as  in  Eq.  4.3  and  y  =  J/ksT.  The 
high-temperature  coefficients,  f{n,j),  are  listed  in  Table  4.9  for  values  of  n  <  8  and  the 
coefficients  /(9, 0), /(9, 1), /(10,0),  and  /(lO,  1)  included  for  completeness.  These  values 
were  calculated  for  a  face-centered-cubic  (fee)  solid;  however,  they  can  be  used  for  an  hep 
solid  because  both  solids  have  12  nearest  neighbors.^ 

The  entropy,  obtained  by  taking  the  temperature  derivative  of  the  free  energy,  is  given 

by 


S^R 


\ni2coshx)-xp+Y,  E  ^-^(^'-^Vfl- 


n>l 0<j<n 


2jx 


pcosh  x 


(4.8) 


The  expression  for  the  pressure  can  be  determined  by  taking  the  derivative  of  the  free 
energy  with  respect  to  the  volume,  giving 


(4.9) 


It  is  easy  to  show  that  the  multiple-exchange  model  with  Jy  =  J2  =  J  and  K'^^  =  0  is 
equivalent  to  the  Heisenberg  model.  For  example,  the  first  order  term  of  Eq.  4.3  can  be 
written  as 


kB 


(4.10) 


^Note  that  the  first  term  in  the  free  energy  is  just  the  Zeeman  energy  term  which  is  also  contained  in 
the  free  energy  for  the  multiple-exchange  model. 

^With  the  assumption  that  the  in-plane  and  out-of-basal  plane  interactions  are  the  same  (Ji  =  J2)  and 
the  four-particle  exchange  process  is  negligible  (A',,  =  0)  in  the  hep  solid. 
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Table  4.9.  High-temperature  expansion  coefficients,  obtained  from  Ref.  [50],  of  the  Heisen- 
berg  model.  For  a  given  n  <  8,  f{n,j)  are  listed  with  increasing  0  <  j  <  ra.  Coefficients  for 
n  =  9, 10  are  not  completely  determined;  however,  available  values  of  /(9, 0), /(9, 1),  and 
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If  one  evaluates  only  the  n  =  1  terms  in  Eq.  4.7,  one  gets 


F  = 


-RTy 


(/(l,0)  +  /(l,iy)  =  -^;,^ 


(4.11) 


which  is  identical  to  that  of  the  exchange  model.  For  n  =  2  the  free  energy  for  the  multiple- 
exchange  model  has  the  form: 


2n  «2(P^) 


-^^(^-^)8(W^ 

^Ji  +  9^2  +  (33J7  +  33J|  +  72Ji J2)P^  '  ■ :  ,    •.  .  8 

■N{l-p^)  ,,,,         ^  =  -A^„,,  ^,(18+120p^  -  138p''). 


(4.12) 
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The  free  energy  of  the  Heisenberg  nearest-neighbor  model  yields 

2 

F  =   -flr|-[/(2,0)  +  /(2,l)pV/(2,2)p*]  (4.13) 

Again,  the  two  models  give  the  same  result  for  a  single  exchange  parameter  and  one  can 
show  that  up  to  n  =  4  the  two  expansions  have  the  same  form  providing  that  Jy  =  J2  =  J 
and  K'g^  =  0  in  the  multiple-exchange  model. 

To  calculate  the  pressure  change  during  a  sample  demagnetization  to  lower  temperature, 
we  employed  the  high-temperature  expansion  to  the  10th  order  in  l/ksT  of  the  Heisenberg 
model.  The  simulation  was  carried  out  in  the  same  manner  as  described  for  the  multiple- 
exchange  model.  The  results,  shown  in  Fig.  4.17,  indicate  a  pressure  increase  at  low  fields  for 
all  entropies.  These  results  validate  our  simulated  data  of  the  multiple-exchange  model  and 
confirm  our  interpretation  that  the  pressure  increase  at  low  fields  of  high-entropy  sample 
is  a  consequence  of  adiabatic  demagnetization  of  a  disordered  ferromagnet. 

High-temperature  expansion  of  the  multiple-exchange  and  Heisenberg  models  can  only 
explain  the  pressure  data  of  high-entropy  samples.  For  both  models,  it  fails  to  predict  the 
pressure  change  of  a  sample  with  low  entropy.  This  would  be  expected  since  the  expansion 
fails  for  T  <  Tc.  We  conclude  that  the  pressure  behavior  observed  for  entropies  of  0.46, 
0.23,  and  0.02i?ln2  are  the  results  of  magnetic  ordering  in  hep  solid  ^He.  In  this  experiment 
the  temperature  of  the  spins  cannot  be  measured  during  the  demagnetization;  however,  the 
ordering  temperature  Tc  can  be  estimated.  Based  on  the  susceptibility  data  [3]  and  the 
calculated  exchange  frequencies  shown  in  Table  4.10  for  hep  solid^,  the  Tc  of  our  sample  is 

*D.  Ceperley,  private  communication. 
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Figure  4.17.  Calculated  pressure  change  using  the  high-temperature  expansion  of  the 
Heisenberg  model  during  a  sample  demagnetization  at  different  initial  sample  entropy. 
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Table  4.10.  Calculated  frequencies  of  exchange  cycles  in  hep  solid  ^He  for  v  =  19.4  cm-'/mole 
by  D.  Ceperley  (private  communication). 


Name 

Frequency  (//K) 

Error  (%) 

J nn 

2.4 

11 

J' 

3.3 

10 

Jt 

4.8 

11 

J 'J' 

2.4 

8 

J ^ 

0.5 

13 

Ka 

1.4 

15 

approximately  17  //K  (assuming  that  Tc  varies  as  v^^^^).  If  the  relation  B/T  =  constant 
for  non-interacting  spins  were  followed,  the  field  at  T  =  17  would  be  2.5(17/614)  = 
0.07  T.  From  Fig  4.8  one  can  see  that  for  B  <  0.07  T  the  pressure  of  the  sample  drops 
rapidly  with  decreasing  field  indicating  a  significant  change  in  the  magnetic  state  of  the 
nuclear  spins.  Thus,  we  conclude  that  Tc  =  17  ±  3  /iK. 


CHAPTER  5 
SUMMARY  AND  FUTURE  WORK 

Theoretical  calculations  [1]  have  predicted  that  the  hep  solid  will  order  ferromagnetically 
with  an  ordering  temperature  of  20  /xK  [23].  Such  a  low  temperature  makes  the  study  of 
the  magnetic  properties  of  the  hep  solid  difficult.  Magnetization  [2,  3]  and  pressure  [4] 
measurements  of  hep  solids  for  T  >  have  indicated  a  strong  ferromagnetic  tendency. 
However,  no  magnetic  transition  was  observed  because  the  ordering  temperature  could  not 
be  reached  using  the  indirect  cooling  method. 

To  achieve  lower  temperature,  ^He  samples  were  directly  demagnetized,  an  idea  that 
has  been  around  for  may  years  [5].  Using  the  sample-demagnetization  technique,  we  carried 
out  pressure  measurements  for  one  bcc  and  two  hep  samples.  With  a  better  design  of  the 
experimental  cell,  we  have  reduced  heating  near  the  end  of  the  sample  demagnetization  and 
eliminated  the  difficulty  encountered  by  the  Tokyo  group.  For  the  first  time,  we  were  able 
to  observe  and  identify  the  signature  of  magnetic  ordering  in  the  hep  solid  -^He.  Below  are 
the  summaries  of  our  pressure  measurements  on  bcc  and  hep  solids.  . .  >    .  < 

5.1    bcc  Solid  •  '  J '  ' 

Our  bcc  sample  with  a  molar  volume  of  20.10  cm^/mole  was  used  as  a  test  of  our  exper- 
imental arrangement.  The  sample  was  precooled  to  0.55  mK  in  a  field  of  2.5  T,  resulting  in 
the  removal  of  nearly  0.99i?ln2  of  the  entropy.  During  the  sample  demagnetization,  we  ob- 
served magnetic  transitions  occurring  when  the  sample  entered  the  HFP  and  the  LFP.  The 
observation  of  the  LFP  transition  indicated  that  our  sample  was  cooled  to  a  temperature 
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below  30  fiK.  Near  zero  field,  the  pressure  of  the  sample  remained  constant,  indicating  a 
much  smaller  heat  flow  into  our  sample  than  in  that  of  the  Tokyo  group  [6]  who  observed  a 
pressure  drop  of  approximately  35  Pa  between  0.01  and  0  T.  This  knowledge  is  important 
in  interpreting  the  results  of  our  hep  samples. 

5.2    hep  Solid 

We  made  pressure  measurements  on  two  hep  samples  with  molar  volumes  of  19.65  and 
19.61  cm^/mole.  There  were  three  demagnetizations  for  the  19.65  cm^/mole  sample  and 
each  was  carried  out  at  different  remaining  entropies  after  a  precooling.  In  the  first  demag- 
netization with  approximately  0.02i?ln2  entropy,  we  found  that  the  pressure  decreased  by 
10  Pa  as  the  field  was  reduced  from  0.2  to  0  T.  Both  the  second  and  third  demagnetizations 
were  with  approximately  0.75i21n2  remaining  entropy  but  at  two  different  temperatures  of 
the  copper  bundle.  In  both  cases,  we  found  a  pressure  increase  in  low  fields  in  contrast  with 
the  0.02i21n2  sample.  For  the  second  demagnetization,  the  pressure  rose  by  26  Pa  between 
0.1  and  0  T. 

The  fourth  demagnetization  was  for  the  19.61  cm^/mole  sample  with  a  remaining  en- 
tropy of  0.85i21n2  and  we  observed  a  pressure  increase  of  20  Pa  between  0.1  and  0  T.  The 
fifth  demagnetization  of  0.46i21n2  entropy  yielded  zero  pressure  change  down  to  zero  field. 
The  sixth  demagnetization  with  0.23i21n2  was  ruined  by  a  large  heat  leak  and  the  seventh, 
a  repeat  of  the  sixth,  yielded  a  3  Pa  pressure  decrease  between  fields  of  0.2  and  0  T. 

Simulated  demagnetizations  of  multiple-exchange  and  Heisenberg  models  indicate  a 
pressure  rise  in  low  field  for  all  entropies.  However,  we  observed  such  an  increase  in  pressure 
only  for  samples  with  0.75/?ln2  entropies  or  higher.  Therefore,  we  conclude  that  the  pressure 
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decrease  observed  in  the  demagnetization  of  a  low-entropy  sample  is  the  ferromagnetic 
ordering  signature  of  the  hep  solid. 

5.3    Future  Work  on  Solid  ^He 

Although  our  pressure  measurements  show  clear  evidence  of  magnetic  ordering  in  hep, 
there  are  other  experiments  which  could  be  done  to  further  our  understanding  of  the  mag- 
netic properties  of  hep  solid.  Two  such  experiments  are  ac  magnetic-susceptibility  mea- 
surement of  hep  solid  and  the  pressure  measurement  of  a  bcc-hcp  mixture,  which  will  not 
only  provide  information  on  the  magnetic  properties  of  the  hep  solid  but  also  of  the  bcc 
solid. ^  .  .  , 

In  addition,  we  suggest  the  measurements  of  Bd  as  a  function  of  v  using  the  direct 
demagnetization  method  for  a  better  determination  of  the  exchange  parameters  of  the  bcc 
solid.  For  the  hep  solid,  we  proposed  the  measurements  of  P„(r,  5o)  at  different  v.  Below 
are  general  description  of  these  experiments. 

5.3.1    AC  Susceptibilitv  Measurement  of  hep  ^He 

The  goal  of  this  experiment  would  be  to  study  the  magnetic  properties  of  hep  solid  in 
its  ordered  state.  Our  pressure  measurements  of  hep  samples  show  that  the  direct  nuclear 
demagnetization  method  can  cool  a  sample  into  its  ordered  state.  This  cooling  method 
therefore  would  be  used  in  the  susceptibility  experiment.  During  the  cooling  and  warming 
processes,  the  magnetic  susceptibility  of  the  sample  would  be  measured  so  that  the  magnetic 
state  of  the  sample  could  be  studied  both  above  and  below  the  transition. 

Although  a  new  experimental  cell  would  be  needed  for  the  susceptibility  measurement, 
the  same  2.5  T  magnet  used  in  the  pressure  measurements  can  be  used  to  provide  the 
'E.  D.  Adams,  P.  Kumeir,  and  Y.  T2ikano,  private  communication. 
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precooling  field  for  the  sample.  The  experimental  cell,  made  from  pure  silver,  would  have 
two  identical  cylindrical  cavities  packed  with  silver  sinter.  One  cavity  would  be  used  as 
the  ^He  cell  and  the  other,  kept  in  vacuum,  would  be  used  as  a  dummy  for  cancelling  out 
the  background  of  the  silver  cell.  The  pick-up  coil,  consisting  of  an  astatic  pair  would  be 
wound  directly  onto  the  body  of  the  cell,  and  the  secondary,  primary,  and  compensation 
coils  would  be  wound  on  a  former  situated  inside  the  static-field  coil. 

In  the  experiment,  an  hep  sample  with  density  near  the  bcc-hcp  line  would  be  formed^ 
and  then  directly  demagnetized  to  the  ordered  state  with  precooling  and  demagnetization 
processes  similar  to  those  described  in  Sections  4.2.1  and  4.2.2.  The  field  dependence  of 
the  magnetic  susceptibility  would  then  be  measured  during  the  cooling  and  warming  of 
the  sample.  The  susceptibility  behaviors  above  and  below  the  ordering  temperature  should 
yield  useful  information  about  the  ordered  state  of  hep  solid  "'He. 

5.3.2    Pressure  Measurement  of  a  bcc-hcp  Mixture 

The  bcc-hcp  phase  line  is  essentially  flat  in  the  P  -T  plane  for  temperatures  near  1  K, 
and  solid  formed  with  a  molar  volume  of  about  19.8  cm'^/mole  would  have  a  mixture  of 
bcc  and  hep  structures  [45].  However,  in  the  vicinity  of  Tjv  and  T^,  the  slope  dP/dT  may 
be  large.  Pressure  measurement  of  the  mixture  would  map  out  the  phase  line,  the  slope  of 
which  would  yield  important  information  on  the  magnetic  characteristics  of  the  mixture. 

If  either  the  bcc  or  the  hep  solid  becomes  magnetically  ordered,  there  will  be  a  large 
change  in  the  pressure  slope  because  of  a  large  reduction  in  the  entropy  of  the  ordered 
^To  have  the  highest  ordering  temperature. 
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Figure  5.1.  Entropy  curves  of  the  bcc-hcp  mixture  for  T/v  >  Tc- 

solid.  The  possible  magnitude  of  this  pressure  change  can  be  estimated  using  the  Clausius- 
Clapeyron  equation 

(5.1) 


dT      Vhcp  -  Vbcc  ' 

where  Skcp  and  Sbcc  are  the  entropy  values  of  the  hep  and  the  bee  solids  respectively  and 
Vhcp  and  Vbcc  are  the  molar  volumes  of  the  solids.  The  difference  Vhcp  -  Vbcc  is  approximately 
-0.1  cm^/mole  at  1  K  [45]  and  decreases  somewhat  as  the  temperature  decreases. 

If  the  bcc  solid  ordered  first  {i.e.  >  Tc),  there  would  be  a  reduction  of  approximately 
0.4i21n2  in  the  entropy  of  the  bcc  solid  [29].  Assuming  that  the  entropy  of  the  unordered 
hep  solid  has  the  same  value  as  that  of  the  bcc  solid  above  T/v,  and  inserting  these  values 
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Figure  5.2.  Entropy  curves  of  the  bcc-hcp  mixture  for  T/v  <  Tc. 
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into  Eq.  5.1,  one  obtains 


dP  0.4Rln2 


-20  Pa//xK. 


(5.2) 


dT      -0.1  cm3/mol 


However,  if  the  hep  solid  orders  first  (Tc  >  Tjv)  then 


dP 
dT 


+20  Pa//iK. 


(5.3) 


Figures  5.1  and  5.2  show  the  entropy  curves  for  the  two  cases.  In  either  case,  the  slope 
change  should  be  close  to  20  Pa//zK  which  is  a  magnitude  larger  than  that  of  the  melting 
curve  near  1  mK.  Such  a  large  pressure  slope  could  be  easily  detected  by  our  current 
experimental  cell. 

The  mixture  could  be  cooled  to  low  temperatures  by  the  direct  demagnetization  method 
similar  to  that  described  in  Chapter  3.  The  pressure  change  versus  time,  P{t),  of  the 
mixture  can  be  monitored  during  and  after  the  demagnetization.  During  the  warm  up 
at  zero  field,  the  P{t)  curve  should  have  a  negative  slope  if  the  bcc  solid  has  a  higher 
ordering  temperature,  and  a  positive  slope  if  the  ordering  temperature  in  the  hep  solid  is 
higher.  Figure  5.3  and  5.4  show  the  predicted  warm-up  pressure  curves  for  both  cases.  The 
ordering  temperature  of  the  mixture  at  zero  field  can  be  estimated  from  known  ordering 
temperatures  at  higher  molar  volumes  using  the  simple  power  law  described  in  Chapter  1. 

5.3.3    Measurement  of      (v)  of  bcc  ^He 

Our  bcc  demagnetization  data  (see  Section  4.1.3)  indicated  a  scaling  factor  of  16.7  for 
Bci{v)  while  results  from  other  works  [12,  6]  yielded  the  values  of  14.5  and  14.7,  respectively. 
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Figure  5.3.  The  predicted  pressure  change  during  a  warm  up  of  the  bcc-hcp  mixture  for 
Tn  >  n. 
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Figure  5.4.  The  predicted  pressure  change  during  a  warm  up  of  the  bcc-hcp  mixture  for 
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This  discrepancy  requires  further  investigation  because  of  its  implication  for  the  volume- 
dependence  of  various  exchange  parameters. 

To  achieve  a  wider  range  of  molar  volumes  than  that  of  the  previous  work  [12],  the 
samples  could  be  cooled  using  the  direct  demagnetization  method.  A  temperature  of  less 
than  30  /xK  can  be  achieved  with  this  cooling  method  allowing  access  to  the  LFP  phase  even 
for  the  highest-density  bcc  solid.  Because  Bd  has  a  T'*  dependence  [13],  most  of  the  sample 
entropy  must  be  removed  before  the  sample  demagnetization  in  order  for  the  sample  to  enter 
the  HFP-LFP  phase  line  at  the  lowest  temperature  where  the  temperature-dependence  of 
Bel  is  small,  see  Fig  4.5.  The  magnetic  state  of  the  sample  could  be  monitored  by  the 
pressure  change  in  the  sample.  A  jump  in  pressure  during  the  demagnetization  as  shown 
in  Fig.  4.4  indicates  that  the  sample  has  crossed  the  HFP-LFP  boundary. 

5.3.4    Measurement  of  Exchange  Parameters  of  hep  "^He 

Because  the  objective  in  this  work  was  to  cool  an  hep  sample  into  its  ordered  state,  we 
did  not  attempt  to  measure  the  volume-dependence  of  the  exchange  parameters  dO^/dv  and 
dK'gg/dv  during  the  sample  precooling  with  high  accuracy.  Although,  our  hep  precooling 
data  are  in  agreement  with  those  of  the  Tokyo  group  [6],  there  is  sizable  scatter  in  the  data 
as  shown  in  Figs.  4.6  and  4.7. 

To  achieve  higher-quality  data,  we  suggest  a  few  modifications  during  the  precooling  of 
the  sample.  First,  the  waiting  time,  at  a  fixed  Tbundle,  should  be  longer  allowing  the  tem- 
perature of  the  sample  to  be  close  to  that  of  the  copper  bundle.  For  example,  in  Table  4.2, 
the  waiting  time  at  Tbundte  =  1-78  mK  was  only  22  hours  causing  a  large  uncertainty  in 
AP;  however,  if  the  waiting  time  were  72  hours  or  so,  the  error  in  AP  would  have  been 
smaller.  Second,  the  excitation  voltage  of  the  capacitance  bridge  should  be  increased  from 
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2.12  Vpp  to  approximately  10  Vpp  allowing  the  sensitivity  of  the  strain  gauge  to  improve 
from  5  Pa  to  2  Pa.  Finally,  the  temperature-dependence  of  the  cell  could  be  reduced  further 
by  having  a  reference  capacitor  made  of  the  same  material  (coin-silver)  and  dimension  of 
the  strain  gauge  and  attached  directly  onto  the  experimental  cell.^ 

5.4  Conclusions 

Direct  nuclear  demagnetization  was  used  to  cool  high-density  solid  ^He  into  the  mag- 
netically ordered  phases.  The  entropy  of  the  samples  was  removed  by  precooling  to  low 
temperatures  in  a  field  of  2.5  T.  The  pressure  changes  in  the  samples  during  precooling 
and  demagnetization  were  detected  by  a  capacitive  strain  gauge.  Two  magnetic  transitions 
were  observed  in  a  bcc  sample  were  observed  with  a  remaining  entropy  of  approximately 
0.01i?ln2.  Pressure  data  of  two  hep  samples  at  various  entropies  indicated  different  behav- 
ior near  zero  field.  Results  from  the  high-temperature  expansion  of  the  Heisenberg  model 
can  explain  only  the  pressure  data  with  remaining  entropies  of  0.72i?ln2  or  higher.  We 
conclude  that  the  pressure  anomalies  observed  in  the  samples  with  entropies  of  0.46,  0.23 
and  0.02i?ln2  were  the  results  of  magnetic  ordering  in  hep  solid  ^He. 

For  further  study  of  magnetism  in  solid  ^He,  we  suggested  two  future  experiments: 
the  ac  susceptibility  measurements  of  a  hep  solid  and  pressure  measurements  of  a  bcc-hcp 
mixture.  Additional  measurements  of  Bci{v)  of  bcc  solid,  and  exchange  parameters  in  hep 
solids  could  also  be  made  using  current  experimental  setup.  These  experiments  would  yield 
valuable  information  about  the  magnetic  properties  of  solid  ^He.  •  • 


Then  both  the  strain  gauge  and  the  reference  capacitor  would  have  the  same  temperature  and  field. 
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